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FRITZ GESZTESY AND RUDI WEIKARD 



Abstract. We provide an overview of elliptic algebro-geometric solutions of 
the KdV and AKNS hierarchies, with special emphasis on Floquet theoretic 
and spectral theoretic methods. Our treatment includes an effective charac- 
terization of all stationary elliptic KdV and AKNS solutions based on a theory 
developed by Hermite and Picard. 



1. Introduction. 

The story of J. Scott Russell chasing a soliton for a mile or two along the 
Edinburgh- Glasgow channel in 1834 has been told many times. It is the start- 
ing point of more than 160 years of an exciting history embracing a variety of deep 
mathematical ideas ranging from applied mathematics to algebraic geometry, Lie 
groups, and differential geometry. In this article we want to tell a certain aspect of 
this story, predominantly from an analytic point of view with special emphasis on 
Floquet theoretic and spectral theoretic methods. Other aspects of this story have 



recently been told by Lax [ 147 1 , Palais [173 , and Terng and Uhlenbeck ]21 



We are interested in the class of evolution equations which permit being cast 
into the form L t = [P, L], where (P,L) is a pair of operators, a so called Lax 
pair, and [P, L] denotes their commutator. Typically, for a fixed operator L, there 
is a sequence of operators P such that L t = [P, L] defines an evolution equation. 
Hence we are actually concerned with hierarchies of evolution equations. Stationary 
solutions of such equations, corresponding to commuting operators P and L, are 
related to algebraic curves and are therefore called algebro-geometric solutions. The 
stationary solutions of higher-order equations in the hierarchy play a decisive role 
in the study of the Cauchy problem for the lower-order equations in the hierarchy. 
Therefore, and because of the connection to algebraic geometry, the stationary 
problem has drawn considerable attention. 

The central object of this article is the class of el liptic algebro-geometric solutions 
of the KdV and AKNS hierarchies (cf. Section [2!^ ) . In a sense, this class represents 
a natural generalization of the class of soliton solutions and enjoys a much richer 
structure due to its connections to algebraic geometry inherent in its construction. 
However, while frequently the algebraic aspects of this construction dominate the 
stage, we will purposely portray a different analytical view often neglected (and to 
some extent forgotten) in this context. Consequently, we will focus in the following 
on the interplay between spectral (and Floquet) theoretic properties of the Lax 
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pairs (cf. Section 2.3) defining the integrable evolution equations in question on the 
one hand, and the construction and properties of the underlying compact Riemann 
surface a ssociated with algebro-geometric solutions on the other. 

In our quest to characterize the class of elliptic algebro-geometric solutions of 
soliton equations in an effective manner, we rely heavily on a marvelous theory 
developed by Fuchs, Halphcn, Hcrmite, Mittag-Leffler, and especially, Picard. We 
have chosen to provide a rather extensive bibliography concerning this classical 
work since some of the cornerstones of this theory seem to have been forgotten, 
and at times appear to be independently rediscovered (cf. Subsection 2.9 and the 
bibliographical remarks at the end of Remark 3.9). In particular, Picard's theorem 
(Theorem 2.7), a key in the aforementioned characterization problem, apparently 
had not been used in the extensive body of literature surrounding elliptic algebro- 
geometric solutions of various hierarchies of soliton equations. Picard's theorem 
suggests one consider the independent variable of the differential equation in ques- 
tion as a complex variable, and study the consequences of assuming the existence 
of a meromorphic fundamental system of solutions. This led to the discovery of 
a connection between the existence of such a meromorphic fundamental system of 
solutions of a linear differential equation Ly — zy for all values of z and algebro- 
geometric properties of L. In particular, it connects the existence of a meromorphic 
fundamental system of solutions of Ly = zy to the integrability of the nonlinear 
equations associated with L via the Lax pair formalism, as we demonstrated in the 
KdV and AKNS cases (cf. Theorems 3.12 and 4.8). As it turns out, Picard was 
describing solutions of differential equations with elliptic coefficients, which repre- 
sent simultaneous Floquet solutions with respect to all fundamental periods of the 
underlying period lattice of the torus in question. Hence Floquet theory, and conse- 
quently, spectral theory, naturally enters when analyzing Lax pairs for completely 
integrable evolution equations and their elliptic algebro-geometric solutions. 

Section H provides an introduction into the KdV hierarchy and its elliptic algebro- 
geometric solutions, reviews the necessary Floquet and spectral theoretic back- 
ground for Hill (Lax) operators, presents Picard's theorem on linear differential 
equations with elliptic coefficients, and especially, provides a historical perspective 
of the subject. Section || is devoted to the KdV hierarchy and its stationary so- 
lutions. The case of rational, periodic, and elliptic stationary KdV solutions is 
described in detail and the role of meromorphic fundamental systems of solutions 
for the corresponding Lax operator associated with stationary KdV solutions is 
underscored. This section culminates in an explicit characterization of all ellip- 
tic, simply periodic, and rational stationary KdV solutions (see J97J, |9£|, [ [234| 
for the original results). Our final Section || describes analogous results for the 
AKNS hierarchy, in particular, it contains an effective characterization of all ellip- 
tic algebro-geometric AKNS solutions (cf. Q for the original proof). We hope the 
enormous bibliography at the end (still necessarily incomplete), will lose some of 
its intimidating character once the reader begins to appreciate the large body of 
knowledge amassed by some of the giants in the field, starting with Hermite. 

Finally we remark that connections between completely integrable sytems and 
Seiberg-Witten theory via elliptic Calogero-Moser-type models, have recently led 
to a strong resurgence of the field of elliptic solutions of soliton equations. It seems 
difficult to keep up with all the current activities in the corresponding preprint 
archives. Hence we refer, for instance, to Babelon and Talon |l^], Donagi and 
Markman Q , Donagi and Witten Q , Itoyama and Morozov ]120| , Krichever ]140[ , 
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Krichever, Babelon, Billey, and Talon [ 141 1 , Krichever and Phong [142], Krichever, 
Wiegmann, and Zabrodin [143], Krichever and Zabrodin [ 144 1 , Kuznetsov, Nijhoff, 
and Sklyanin 145 1 , Levin and Olshanetsky [ 149 1 , Marshakov [152], and Vaninsky 
229 ), from which the interested reader can easily find further sources. 



2. Background and some History. 

2.1. The Korteweg-de Vries Equation. In the first few decades after their 
discovery, solitary waves were considered, for instance, by Stokes, Boussinesq, and 
Lord Rayleigh. But the most far reaching, and in a sense, lasting contribution 
was made in 1895 by Korteweg and de Vries [129], who deduced their celebrated 
equation, which may be cast in the form 



(2.1.1) 



1 



-Mx 



(However, as recently alluded to by P ego [178], Boussinesq originally noted a sys- 
tem of equations equivalent to ( 2.1.1 ) and used it to study solitary waves in the 
1870's.) In 1960 Gardner and Morikawa j79| used the KdV equation to describe 
collisionless-plasma magnetohydromagnetic waves. Since then the KdV equation 
has been rederived again and again in different contexts as a model equation de- 
scribing a considerable variety of physical phenomena, and can now be considered 
one of the basic equations in mathematical physics. Even more important is the 
fact that this development triggered the examination of a whole class of other non- 
linear evolution equations, some of which are of great physical relevance, such as 
the nonlinear Schrodinger equation, the sine-Gordon equation, the Toda lattice, the 
Boussinesq equation, the Kadomtsev-Petviashvili equation, etc. 

But the description of wave phenomena alone would not have turned the subject 
into the kind of industry it is today. In 1955 Fermi, Pasta, and Ulam [[37] studied 
a system of nonlinear oscillators which may be viewed as a discretized version 
of the KdV equation on a finite interval with periodic boundary conditions. To 
their surprise they found that energy shows little tendency toward equipartition 
among the degrees of freedom. Later Kruskal and Zabusky 1 241 examined the 
KdV equation numerically and observed the formation of solitary waves which, 
using their words, "pass through one another without losing their identity" . In 
order to emphasize this particle-like behavior Kruskal and Zabusky coined the term 
soliton. 



In 1968 Miura [ 164 1 introduced the transformation 

(2.1.2) q(x, t) = z — v(z, x, t) 2 — v x (z, x, t) 

in which z is a (generally complex) spectral parameter. This transformation, now 
known as Miura's transformation, relates the KdV equation to a variant of the so 
called modified KdV (mKdV) equation, 

(2.1.3) v t = v xxx - 6v 2 v x + 6zv x . 



Since (2.1.2) is a Riccati equation for v, it may be transformed into the linear 
equation 

(L(t)y)(z, x, t) = y"(z, x, t) + q(x, t)y(z, x, t) = zy(z, x, t) 

by introducing v(z,x, f) = y'{z, x, t)/y(z, x, t) (primes denote derivatives with re- 
spect to x). Gardner, Kruskal, and Miura Jl65| showed that the eigenvalues of the 
£ 2 (R)-operator associated with L{t) do not depend on t, that is, they are constants 
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of the motion under the KdV flow. This was the starting point for another seminal 
work, this time by Gardner, Greene, Kruskal, and Miura ]77t , in which they used 
the inverse scattering method to solve the Cauchy problem for the KdV equation 
with rapidly decaying initial data. This method, which represents a nonlinear ana- 
log of the Fourier transform to solve linear partial differential equations, initially 
consists of computing the scattering data of L(0), then propagating them in time 
(which is simple), and finally reconstructing the potential q{x,t) in L(t) via the 
Gelfand-Levitan, or rather, the Marchenko equation. Miraculously, this function 
q(x, t) is the desired solution of the Cauchy problem. See, for instance, Ablowitz 
and Clarkson JD, Ch. 2, Asano and Kato [fLl|, Chs. 5,6, Dodd, Eilbeck, Gibbon, 



and Morris 46|, Ch. 4, Drazin and Johnson ]52], Ch. 4, Gardner, Greene, Kruskal 



and Miura |7q|, Iliev, Khristov, and Kirchev 117, Ch. 3, Lax [ 147 1 , Marchenko 



153 , Ch. 4, and Palais |173 for more details 



2.2. Hamiltonian Systems. The dynamics of a classical mechanical particle sys- 
tem is described as a flow on a symplectic space, called the phase space. The phase 
space is the cotangent bundle of a Riemannian manifold, called the configuration 
space, describing the positions of the particles. The flow is given by Hamilton's 
equation q t — grad s (-ff), where gradg is a symplectic gradient and H is the so 
called Hamiltonian of the system. If there are only finitely many, say n, degrees 
of freedom (as in a system of finitely many interacting particles) then q — (£,??), 
where £ is a vector of local coordinates on the configuration space (called general- 
ized coordinates) and rj is the vector of the associated momenta (called generalized 
momenta). H is then a function of £, ?y, and t, and the symplectic gradient is given 

by 

grad s = ^ ^ grad ( £ !7)) , 

where /„ is the n x n identity matrix. A canonical transformation is a change of 
coordinates in phase space which leaves the form of Hamilton's equation invari- 
ant. Sometimes there exists a canonical transformation such that the transformed 
Hamiltonian is a function of the transformed generalized momenta alone which 
are then called action variables while the transformed generalized coordinates are 
called angle variables. The action variables are then constants of the motion while 
the angle variables change linearly with time. If this happens the system is called 
completely integrable. 



In 1968 Gardner, Kruskal, and Miura [16!:] explicitly constructed an infinite 



sequence of constants of the motion providing the first hint toward a Hamiltonian 



structure of the KdV equation. Next, in 1971, Faddeev and Zakharov [ 242 1 , in 
an attempt to explain the unusual behavior of the KdV equation, showed that it 
can be viewed as an infinite-dimensional completely integrable Hamiltonian system. 
Take, for instance, the Schwartz space <S(R) as phase space which can be viewed as 
a symplectic manifold. A symplectic gradient on S(R) is given by 

grad s = grad , 

where grad(-F) denotes the Lagrangian (or variational) derivative of F 
More precisely, if F is of the form 



/oo 
dxF(q, </,... ,<?(">), 
-OO 
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where F : R n+1 — > K is a polynomial function without a constant term, then 

dF fdF\' ( dF x {n) 

+ - + <- 1 > "Uw 



3FY f dF\" . ._/ 



and hence 

grad g (F)=^j +... + (-1) 

using primes to denote derivatives with respect to x. The Hamiltonian flow 

q = grad s (F) 

then becomes the nonlinear evolution equation 

fdF\' fdF\" , / dF \ {n+1) 

In the special case of the KdV flow, the Hamiltonian H is given by 

H:S(R)^R, q^ij dx(q(x) 3 --q'(x) 2 ), 

and Hamilton's equation qt — gr&d s (H) results in the KdV equation ( |2.1.l| ). Anal- 
ogous considerations apply to the periodic case replacing dx by J Q dx, with 
fl > the fundamental period of q. We refer, for instance, to Palais |173| for more 
information. 

2.3. Lax Pairs and the KdV Hierarchy. Gardner, Greene, Kruskal, and Miura 
77 had shown that the L 2 -spectrum of Lit) — d 2 /dx 2 + q(x,t) is independent of 



t whenever q is a solution of the KdV equation. This discovery inspired Lax to 
conjecture that the operators L(t) are all unitarily equivalent to one another, that 
is, L(0) = U (t)^ 1 L(t)U (t) for some family of unitary operators U(t). This then led 



to his celebrated commutator representation [146], 

q t = L t = [P 3) L], 

where P3 denotes the third-order ordinary differential expression defined by 

D _ 1 d 3 3d 3 

which governs the time evolution of U(t) according to U t (t) — Ps(t)U(t), U(Q) = I. 
Indeed, the latter equation implies (U(t)~ 1 ) t = —U(t)~ 1 P 3 (t) and hence formally, 

d 



yields 



dt (U(t)-'L(t)U(t)) = U{t)-\L t {t) - [P 3 (t),L(t)])U{t) = 



Uty-^LtyUit) = L(0), t G 



(Functional analytic arguments guaranteeing self-adjointness of L(t) and unitarity 
of U{t) in the Hilbert space L 2 (R;dx) can easily be supplied.) The letters P and 
L were first used by Gelfand and Dickey [|4| in honor of Peter Lax. Accordingly, 
(P3, L) is called a Lax pair. 

In fact, Lax designed his procedure for a far more general setting: it applies to 
evolution equations for q whenever there is a one-to-one correspondence between 
q(-,t) and a (formally) symmetric operator Lit), The goal is then to find a (for- 
mally) skew-symmetric operator P(t) satisfying U t {t) = P(t)U(t), where L(0) = 
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U [t)^ 1 L(t)U (t) . In particular, returning to the case where Lit) = d 2 jdx 2 + q(x, t), 
Lax showed that one can always find an odd-order differential expression P 2n +i(t) 
such that [P 2n +i(t), L(t)] is an operator of multiplication. The coefficients of 
P2n+i(t) are then differential polynomials in q(x,t), that is, polynomials of q(x,t) 
and its rr-derivatives. The commutator [P 2n +i, L] is then also a differential polyno- 
mial in q and the evolution equation 

(2.3.1) (? t = [P 2n+1) L] 

defines the n-th equation of the KdV hierarchy. 

More precisely, if L = d 2 /dx 2 + q, then any ordinary differential expression P, 
for which [P, L] equals an operator of multiplication, is the sum of a polynomial 
K € C[L] and an odd-order differential expression P2„+i, that is, 

P = K(L) + P 2n+1 . 

Upon rescaling the time variable, P 2n +i can be assumed to be monic. Moreover, 
choosing the polynomial K appropriately, it may be shown that 



(2.3.2) P 2n+1 = ^ 

j=0 



--f'+f- 



for some integer n e No(= NU {0}), where the functions fj satisfy the recursion 
relation 

fo = 1, 

(2.3.3) /;., [c ■ <//; • ./ " "■ 

Then 

[P, L] = [P 2n+1> L] = -/;" + qf' n + ±q'f n 

and introducing 



P n (z,x,t) — / J fn—j (•£? t)z^ , 
3=0 

this recursion relation becomes 

(2.3.4) [P 2n +i,L] = ip/" + 2(q - z)F n ' + q'F n . 

Conversely, if P„ (z, x) is a polynomial in z such that F n "' + A(q — z)F^ + 2q'F n does 
not depend on z, then its coefficients define a differential expression P 2n +i such 
that equation (2.3.4) is satisfied. In particular, if F n '" + 4(g — z)F' n + 2q'F n = 0, 
then q is a stationary (i.e., t-independent) solution of one of the n-th equations in 
the KdV hierarchy. 

The first few equations of the KdV hierarchy explicitly read, 

Qt = q x , 

1 3 

qt = -^q_xxx + -^qq x + ciq x , 

1 5 5 15 2 ,\ 3 , 

qt = -^qxxxxx + ^qqxxx + -^q x qxx + -g-g q x + ci{-q xxx + -qq x ) + c 2 q x , 

etc., 
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where the ci denote integration constants encountered in solving the recursion 
( 2.3.3| ) for fj+i- Here we ought to mention that each of these equations is a com- 
pletely integrable Hamiltonian system and that the sequence of these equations 
is intimately related to the sequence of conservation laws discovered by Gardner, 



Kruskal, and Miura [165]. 

However, the Lax method can be further generalized in a variety of ways: If 
one considers instead of I a first-order 2 x 2-matrix differential expression (i.e., a 



Dirac-type operator), one arrives at the ZS hierarchy [243 which is associated with 
the nonlinear Schrodinger equation, and more generally, at the AKNS hierarchy Q 
which includes the KdV, the mKdV, and the nonlinear Schrodinger hierarchies as 
special cases. If one considers instead of the second-order differential expression 
L a differential expressions of any order, one arrives at the Gelfand-Dickey hierar- 
chy ^5|, j84p. Moreover, the formalism is not confined to differential expressions 
but extends to formal pseudo-differential expressions (including rational functions 
of d/dx) needed in connection with the sine-Gordon and (modified) Kadomtsev- 
Petviashvili hierarchies. 

2.4. The Algebra of Commuting Differential Expressions. As seen from 



( 2.3.1), stationary (i.e., time-independent) Lax equations naturally lead to commut- 
ing differential expressions. Independently of this fact, the question of commuting 
differential expressions was raised by Floquet J7(J in 1879. Some 25 years later, 

To this 



the question was again considered by Wallenberg p31| and Schur |19J 
end, Schur developed the algebra of symbolic, that is, formal pseudo-differential 
expressions. He proves, in particular, the following statement: if A and B are two 
monic commuting differential expressions then the coefficients of A are differential 
polynomials in the coefficients of B (for a contemporary approach to these results, 
see, e.g., Wilson [ 236 1 ) . The decisive step, however, was done by Burchnall and 
Chaundy in the 1920's. They proved the following result in p0[ . 

Theorem 2.1. Let A and B be ordinary differential expressions of relatively prime 
orders m and n, respectively. Then [A, B] = if and only if there exists a polynomial 
f of the form 



(2.4.1) f(a,0) =a n -{3 m + £ c j>k oj p k , 

such that f{A,B) = 0. 



j,k>Q 
mj+nk<mn 



f is called the Burchnall-Chaundy polynomial of A and B. In |31J] Burchnall 
and Chaundy constructed differential expressions A and B satisfying the equation 
f(A,B) = for a given polynomial / of the form ( 2.4.1 ) (cf. also Baker jl6|, 
Burchnall and Chaundy |52|). More recent treatments of the Burchnall-Chaundy 
theory can be found, for instance, in Carlson and Goodearl |37|1, Gatto and Greco 



134), Mumford |6§, Previato |86|, and Wilson J237 ] 



, Giertz, Kw ong , and Zettl J 102| , Greco and Pre viato [104 , Krichever [ 133 1 , 



Finally, we briefly return to the Lax pairs (P,L) discussed in Section 2.3. Be- 
cause of the Burchnall-Chaundy relationship we call L (or the set of its coefficients) 
algebro-geometric if there exists a corresponding P such that [P, L] = 0. We will 
provide more precise definitions in the contexts of the KdV and the AKNS hierar- 
chies in Definitions |3 . l| and |4.2|, respectively. 
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2.5. Elliptic Functions in a Nutshell. Elliptic functions provide some of the 
most important examples of algebro-geometric potentials. We present here a very 
brief account of Weierstrass' point of view. For general references see, for instance, 



Akhiezer M, Chandrasekharan |3§J], Markushevich [154], and Whittaker and Wat 



son 1 235] 



A function / : C — * C U {oo} with two periods a and 6, the ratio of which is not 
real, is called doubly periodic. If all its periods are of the form m\a + TO2& where 
mi and m,2 are integers then a and b are called fundamental periods of /. 

A doubly periodic meromorphic function is called elliptic. 

It is customary to denote the fundamental periods of an elliptic function by 2u>\ 
and 2^3 with Im(o>3/o>i) > 0. We also introduce ui2 — w\ + u>3 and o>4 — 0. The 
numbers o>i, ...,o>4 are called half-periods. The fundamental period parallelogram 
(f.p.p.) A is the half-open region consisting of the line segments [0, 2o>i), [0,20*3) 
and the interior of the parallelogram with vertices 0, 2o>i, 2o>2 and 2o>3. 

The class of elliptic functions with fundamental periods 2o>i , 2o> 3 is closed under 
addition, subtraction, multiplication, division by non-zero divisors and differentia- 
tion. If / is an entire elliptic function then / is constant. A non-constant elliptic 
function / must have at least one pole in A and the total number of poles in A 
is finite. The total number of poles of an elliptic function / in A (counting mul- 
tiplicities) is called the order of /. The sum of residues of an elliptic function / 
at all its poles in A equals zero. In particular, the order of a non-constant elliptic 
function / is at least 2. The total number of points in A where the non-constant 
elliptic function / assumes the value A (counting multiplicities), denoted by n(A), 
is equal to the order of /. In particular, n(A) > 2. Furthermore, s(A), the sum of 
all the points in A where the non-constant elliptic function / assumes the value A, 
is congruent to s(oo), the sum of all the points in A where / has a pole, that is, 
s(A) — s(oo) + 2mio>i + 2TO3O3, where m\ and 7713 are certain integers. 

The function 

p(z;u> 1 ,U 3 ) = ^ + ( ( z - 2mwi - 2no> 3 ) 2 ~ (2mwi + 2no> 3 ) 2 

(m,»')^(0,D) 

or p(z) for short, was introduced by Weierstrass. It is an even elliptic function of 
order 2 with fundamental periods 2o>i and 2o>3. Its derivative p' is an odd elliptic 
function of order 3 with fundamental periods 2o>i and 203. Every elliptic function 
may be written as R±(p{z))+R2{p(z))p'{z) where R\ and R2 are rational functions 
of p. 

The numbers 

92 = 60 ^ ' 



m,ne: 

(m,n)#(0,0) 



.93 = 140 Y, 



(2mo>i + 2tmj 3 ) 4 ' 
1 



(2mo>i + 2no> 3 ) 

(m,nj#(0,0) 



are called the invariants of p. Since the coefficients of the Laurent expansions of 
p(z) and p'(z) at z — are polynomials of 172 and 33 with rational coefficients, the 
function p(z; o>i, 0/3) is also uniquely characterized by its invariants 172 and 53. One 
frequently also uses the notation p{z\g2,g%). 



ELLIPTIC SOLUTIONS OF SOLITON HIERARCHIES 



9 



The function p(z) satisfies the first order differential equation 

(2.5.1) p'(z) 2 =4p(z) 3 -.g 2 p(z)-« ?3 

and hence the equations 

p"{z) = 6p(z) 2 - g 2 /2 and p'"{z) = I2p\z)p{z) 

which shows that — 2p is a stationary solution of the KdV equation. 

The function p', being of order 3, has three zeros in A. Since p' is odd and 
elliptic it is obvious that these zeros are the half-periods uj\ , w 2 — ^i + ^3 and 
ui 3 . Let €j = p(u>j), j = 1, 2, 3. Then ( [2.5.1 ) implies that 4e? — g 2 ej — g 3 = for 
j = 1,2,3. Therefore 

= ei + e 2 + e 3 , 

92 = -4(eie 2 + eie 3 + e 2 e 3 ) = 2{e\ + e\ + e 3 ), 
.g 3 = 4eie 2 e 3 = ^(e? + e 2 ! + e 3 ; ). 

Weierstrass also introduced two other functions denoted by C and a. The Weier- 
strass ^-function is defined by 

ICW = -**), iim(C(,)-i) = o. 

It is a meromorphic function with simple poles at 2uiloi + 2nw 3 , m, n G Z having 
residues 1. It is not periodic but quasi-periodic in the sense that 

C(* + 2ljj) = C(z) + 2r)j, j= 1,2,3,4, 

where rjj — ({loj) for j = 1,2,3 and 774 = 0. 
The Weierstrass cr-function is defined by 

— = COO, lim — = 1. 

a is an entire function with simple zeros at the points 2mLO\ + 2nw^,m,n G Z. 
Under translation by a period a behaves according to 

a[z + 2u j ) = -a{z)e 2 ^ z+w ^ , j = 1,2,3. 

Next we recall the following fundamental theorems. 

Theorem 2.2. Given an elliptic function f with fundamental periods 2uj\ and2u>3, 
let 61, fe r 6e i/ie distinct poles of f in A. Suppose the principal part of the Laurent 
expansion near bk is given by 



0H , 

3= 



k = 1, ...,r. 



/(*) = c + E E(- 1 )^ 1 7^tvC ( - 1) (- - *o, 

where C is a suitable constant and( is constructed from the fundamental periods 2u>i 
and 2u>3. Conversely, every such function is an elliptic function ifY^k=i -^1,* = 0- 
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Theorem 2.3. Given an elliptic function f of order n with fundamental periods 
2lj\ and 2loj,, let ai,...,a n and b±,...,b n be the zeros and poles of f in A repeated 
according to their multiplicities. Then 

v _ a(z - ai) ■ ■ ■ a(z - a n ) 

I[Z) ~ a{z - b x ) ■ ■ ■ a(z - K-Mz - b' n ) ' 
where C is a suitable constant, a is constructed from the fundamental periods 2lo\ 
and 2lo3 and where 

K -b„= ( a i + ■•■ + a n) - (h + ■■■ + b n ) 

is a period of f. Conversely, every such function is an elliptic function. 

Finally, we turn to elliptic functions of the second kind, the central object in 
our analysis. A meromorphic function ip : C — > C U {00} for which there exist two 
complex constants toi and u>3 with non-real ratio and two complex constants pi and 
P3 such that for i = 1,3 

ip(z + 2ui) = Piip(z) 

is called elliptic of the second kind. We call 2u>i and 2oj 3 the quasi-periods of ip. 
Together with 2u>i and 2w 3 , 2m 1 u; 1 + 2m 3 w 3 are also quasi-periods of ip if mi and 
7713 are integers. If every quasi-period of ip can be written as an integer linear 
combination of 2u>i and 2^3 then these are called fundamental quasi-periods. 

Theorem 2.4. A function ip which is elliptic of the second kind and has funda- 
mental quasi-periods 2u)\ and 2lut, can always be put in the form 

if \ n r \ ,cr(z - a x ) ■ ■ ■ a(z - a n ) 

ip(z) = Cexp(Az) — — - — 

cr(z - bi) ■ ■ ■ a{z - b n ) 

for suitable constants C, X, a\, ...,a n and b\, ...,b n . Here a is constructed from the 
fundamental periods 2ui\ and 2u>3 . Conversely, every such function is elliptic of the 
second kind. 

Theorem 2.5. Given numbers a\, ...,a m and (3\, m such that (3k 7^ Pe (mod A) 
for k £, the following identity holds 

-pj tr(x-aj) = y> rifcLi o{fh - «fc) (t{x ~ (i 3 + (3 - a) 
1 ' ■ ] 1 = 1 a(x ~ fa) £j Uti,i^ " ft) ~ 0iM0 ~ «) ' 

where 

m m 

a = otj and (3 = (3j 
3=1 j=i 
and a is constructed from the fundamental periods 2oj\ and 2lo 3 . 

Sketch of proof. Since this result seems less familiar than those above we will briefly 



sketch its proof. Denote the left and right-hand sides of (2.5.2) by / and g, respec- 
tively. Both / and g are elliptic functions of the second kind associated with the 
same Floquet multipliers (with respect to translations by 2tUk, k = 1, 3). Their quo- 
tient is therefore an elliptic function. Also / and g have the same poles and zeros 
taking multiplicities into account. The statement about the zeros is a consequence 
of the identity (cf. (235], p. 451) 



in / T~f 7n ( Si \ \ m " L 

<«■»> 
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Therefore f/g is entire and hence constant. Since all residues are equal to one 

/ = <?• □ 

2.6. Hill's Equation and its Spectral Theory. The study of linear homoge- 
neous differential equations with periodic coefficients predates the late nineteenth 
century, but the equation (Ly)(x) = y"(x) + q(x)y{x) = zy(x), where q(x) is a 
continuous, real- valued, periodic function of a real variable x and z is a real param- 
eter, has generally been called Hill's equation since its appearance in the study of 



the lunar perigee by Hill [114| in 1877. In addition to its applications in celestial 
mechanics, this equation has found countless applications in quantum mechanics, 
where it becomes Schrodinger's equation and is used, for instance, to model crystal 
structures of solids. 

Periodic differential equations are usually studied by applying Floquet theory. 
Floquet theory (first developed by Floquet ]7l] |, J72[ starting in 1880) specifies the 
general structure of solutions of systems of periodic differential equations. Consider 
the equation y'(x) — Q(x)y(x), where Q is an n X n matrix whose entries are 
continuous (for simplicity) and periodic with period O > and y{x) is C"-valued. 
Let y be the space of solutions of y'{x) = Q(x)y(x) and Tq the restriction of 
V l— * ?/(■ + to y. Floquet theory then amounts to the study of the operator Tq. 
Since 7q maps the n-dimensional vector space y to itself, the problem is reduced 
to a problem in linear algebra. The eigenvalues and eigenfunctions of are called 
Floquet multipliers and Floquet functions, respectively. For general references on 
Floquet theory see, for instance, Arscott pOj, Coddington and Levinson p3|, Ch. 3, 



Eastham (5§, Ince pS] , Sect. 10.8, Magnus and Winkler pl| , Marchenko [ 153 1 , 
Sect. 3.4, McKean and van Moerbeke |160| , and Yakubovich and Starzhinskii [ 240 1 . 
In the special case of Hill's equation (Ly)(x) = y"{x) + q{x)y(x) = zy(x), we will 
denote the translation operator restricted to the set of solutions (y(z 7 x),y'(z, x))' e 
y(z) of the associated first-order system, by Tq(z). 

Several differential operators (resp., boundary value problems) are studied in 
connection with Hill's equation, 

(Ly)(x) = y"(x) + q{x)y(x) = zy(x), q{x + Q) = q(x). 

(1) The maximally defined operator H in L 2 (R) associated with L. 

(2) Auxiliary operators in L 2 ([xq, xq + associated with L and certain families 
of boundary conditions, in particular, Dirichlet boundary conditions. 

(3) The operators in L 2 ([xo, xq + ^]) associated with L and cyclic boundary con- 
ditions y(xo +0) = exp(i6)y(x ) 7 y'(xo +0) = exp(i9)y' (xo) , 9 G [0, 2n), in 
particular, periodic (0 = 0) and anti-periodic (9 = tt) boundary conditions. 

Floquet theory provides us with a handle on the problem of determining spectral 
properties of these operators: the periodic eigenvalues are given as the (necessarily 
real) zeros of tr(7o(z)) — 2, while the anti-periodic eigenvalues are given as the 
(necessarily real) zeros of tr(Tn(z)) + 2, and the spectrum of H coincides with 
the conditional stability set S(q), that is, the set of all values z € K such that 
(Ly)(x) = zy(x) has a nontrivial bounded solution with respect to x S K. This was 
first shown by Wintner ]238| , [ |239[ in 1947/48. The conditional stability set, in 
turn, may be characterized as the set of all z £ R such that Ly — zy has a Floquet 
multiplier of absolute value one. Since 



det(T n (z)) = 1, 
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the Floquet multipliers, being the eigenvalues of 7q(z), are given as the zeros of 
p 2 — /otr(7h(z)) + 1 and we get 

S(q) = {z G R | - 2 < tr(T (z)) < 2}. 

The conditional stability set consists of countably (possibly finitely) many closed in- 
tervals (plus possibly a half- line) , whose endpoints coincide with points where only 
one (linearly independent) Floquet solution exists (which is necessarily (anti-)perio- 
dic). This was first shown by Hamel JllO[ in 1913 (see also Liapunov ]150| , who 
treated the case y" = Xpy with a periodic function p in 1899, and Haupt |lllj , who 
corrected a mistake in Hamel's paper). Hence the (anti-)periodic eigenvalues deter- 
mine the spectrum of H. From the periodic and anti-periodic eigenvalues repeated 
according to their multiplicity and ordered as a decreasing sequence denoted by 
i?o,-E<i, ■■• (observing that for normal operators the algebraic and geometric multi- 
plicities of eigenvalues coincide), one infers 

oc 

a(H)=[j[E 2j+1 ,E 2j }. 

3=0 

Typically, the spectral bands [E^+i, E 2 j] are separated by spectral gaps {E 2 ,j +2l 
i?2j+i)- However, since the (anti-)periodic eigenvalues may be twofold degenerate, 
some gaps (or even all gaps) may close. If only finitely many gaps are present, one 
calls q, the potential coefficient in L, a finite-gap, or a finite-band potential. 

A trivial example of a finite-band potential is of course the constant potential 
q(x) = c. The first in this century to discuss a nontrivial example was Ince [ p_18| 
around 1940, who treated in depth Lame's potential 

(2.6.1) q(x)=-s(s + l)p(x + u 3 ), s G N, 

with fundamental half periods u>\ G R and L03 £ iR. Under these conditions q 
is real-valued and real-analytic in x. Ince showed that for all but 2n + 1 values 
of z G R, the equation y"(x) + q(x)y(x) = zy(x) has two linearly independent 



(anti-)periodic eigenvalues, that is, q in (2.6.1) is a finite-band potential (see also 



Akhiezer ||, Erdelyi |36j, Turbine r [f228[ , and Ward [232]). However, a closer look 
at the classical works of Hermite ]113[ and Halphen [10E] (see also Klein ]125[ ) on 
Lame's equation at the end of last century reveals that these results were actually 
well-known (but not yet put in a Floquet-theoretic language) as can be inferred, 
for instance, from Whittaker- Watson's treatise ]235| ], Sect. 23.7. 

In 1909 Birkhoff |23|] compared eigenvalues of various boundary value problems 
associated with L on a finite interval. His results show that there is precisely one 
Dirichlet eigenvalue in each of the intervals [E 2 j+ 2 , E 2 j+i] (i.e., in the closure of the 
gaps). In contrast to the (anti-)periodic eigenvalues, the Dirichlet eigenvalues asso- 
ciated with the interval [xo , xq +0] vary with xq in [£'2^+2 , -E^j+i] if i?2j+2 < ^2j+i • 
Writing M = U„ e z[^o + nCl, xq + (n + 1)0], a comparison of the Dirichlet problems 
on [xq,Xq + CI) and on (—co,x a ) U (a;o:°o) assuming q(x) to be O-periodic, yields 
the same Dirichlet spectra in the spectral gaps (E 2 j+2, E 2 j+i) with E 2 j+ 2 < ^2j+i- 
Explicit formulas for the monotone rate of change of various kinds of eigenvalues 
(including Dirichlet eigenvalues) with respect to varying xq G R can be found in 



Kong and Zettl [ 126 1 , [[127) and the literature therein. In |89j this phenomenon is 
related to Green's functions and rank-one perturbations of resolvents. 

If q is a locally integrable, real-valued, periodic function with period £1 > we 
therefore have equivalence of the following statements: 
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(A) q is a finite-band potential. 

(B) For only finitely many values of z the differential equation y"(x) +q(x)y(x) = 
zy(x) fails to have two linearly independent Floquet solutions. 

(C) The Dirichlet boundary value problem on the interval [xq, xq + CI] has only a 
finite number of eigenvalues depending on Xq. 

Spectral theoretic aspects of c omplex-valued pe riod ic potentials q have been 
investigated by Rofe-Beketov [18S] and Tkachenko [ 21 1 1 (see also Birnir pij , p5[ , 
Kotani pi, McGarvey p57[-|l59|, Sansuc and Tkachenko [l90|-]T9^|, and Tka- 
chenko [ 2 1 2 1 , [213], 1 215] for recent results). They found that the spectrum of H 
and the conditional stability set of L still coincide, that is, 

S(q) = {z G C | - 2 < tr(T n (z)) < 2}. 

Since the conditional stability set is given as the preimage of [—2, 2] under an entire 
function, it turns out that the spectrum of H consists of countably many (possibly 
finitely many) regular analytic arcs. While the term "finite-gap potential" is now 
rendered meaningless, it still makes perfect sense to call a potential a finite-band 
potential if the spectrum is a finite union of regular analytic arcs. In contrast to 
the real-valued case, these spectral arcs can now cross each other, see, for instance, 
(93), [176 1 for explicit examples exhibiting this phenomenon. 

Returning to the real- valued case, the isospectral set I(qo) of a given periodic 
potential qo G C(R) with period Cl > (i.e., the set of all Sl-periodic q G C(R) whose 
2f2-periodic eigenvalues coincide with that of qo) turns out to be a manifold, in fact, 
a (generally infinite dimensional) torus generated by a product of circles. Each circle 
is uniquely associated with a spectral gap (£2^+2, £2j+i), £2.7+2 < £2,7+1 and the 
periodic motion of the Dirichlet eigenvalue in the gap (£?2j+2, £2,3+1) corresponding 
to the interval [x, x + O] as a function of x G R. In the special case where q is real- 
valued and periodic, and the £ 2 -spectrum of —d 2 /dx 2 + q(x) is a half-line (—00, Eq], 
Borg [^J proved the celebrated inverse spectral result q(x) — Eq for all x G R. A 
quick proof of this uniqueness result and extensions to real-valued reflectionless 
potentials q with associated half-line spectrum (— 00, Eq] follows from the trace 
formula proved in ]89| (observing £(A, x) — 1/2 for A < Eq in the corresponding 
trace formula (3.1) for q(x) in J89|). Similarly, if q is real-valued and periodic, and 
the £ 2 -spec trum of —d 2 /dx 2 + q(x) is of the type (—00, £2] U [£i,£o], £ 2 < E\, 
Hochs tadt [ 116 proved that q(x) — —2p(x + U3 + a), that is, the Lame potential 
( [2.6.1 ) for n = 1, where p(x) denotes the Weierstrass elliptic function associated 
with some period lattice u>i > 0, —iu>3 > (depending on Eq,E\, and E2) and 
a G R. The isospectral torus for real-valued periodic potentials with three or more 
spectral bands is described, for instance, in Buys and Finkel [|i|, Finkel, Isaacson, 
and Trubowitz ]6q|, Gesztesy, Simon, and Teschl M], Gesztesy and Weikard fl92fl, 



Iwasaki [124], McKean van Moerbeke [ 160 1 , and McKean and Trubowitz [161]. 

It should be emphasized in this context that the assumption of real-valuedness 
of q{x) cannot be dropped as shown by the well-known example q(x) — e xp(z x), 
with associated L 2 -spectrum (—00, 0] (see the paragraph following Remark 3.5 for 
more details). 



2.7. Periodic KdV Potentials. In 1974, Novikov |169| ] investigated the Cauchy 
problem of the KdV equation in the case of periodic initial data. He noted that 
the right generalization of multi-soliton solutions, which are stationary solutions of 



14 



FRITZ GESZTESY AND RUDI WEIKARD 



appropriate higher-order KdV equations, are the finite-band potentials, that is, he 
proved the following result. 

Theorem 2.6. If q is a real-valued, periodic, stationary solution of an n-th order 
KdV equation, then the L 2 (M) -spectrum associated with d 2 /dx 2 + q(x) has at most 
n finite bands. 

Within a year, Dubrovin |33) and Flaschka j6^[ also proved the converse, and 
Dubrovin and Novikov (57) used these results to solve the Cauchy problem of the 
KdV equation in the case of finite-band initial data. Hence the statement, 

(D) q is a stationary solution of a higher-order KdV equation, 
is equivalent to any of the statements (A)-(C) in Section 2.6 for real-valued q. 

Moreover, it became clear from these investigations that generically, finite-band 
potentials will only be quasi-periodic with respect to x and not periodic in x. 

About the same time, Its and Matveev [ |123| derived their celebrated formula for 
finite-gap solutions q(x,t) in terms of the Riemann theta function associated with 
the underlying hyperelliptic curve and a fixed homology basis on it (cf. |l7j] , Ch. 3, 
|87| , Ch. 1, [ |170| , Ch. II). Subsequent extensions of this formula to general (matrix- 
valued) integrable systems were developed by Dubrovin f34||-|5l| and Krichever 
|133| . [ |134| , j 137 1 . A new approach to finite-band solutions of the KdV hierarchy in 
terms of Klcinian functions was recently developed by Buchstaber, Enol'skii, and 
Leykin §§, j§. 

2.8. Elliptic KdV Potentials. While the considerations of the preceding subsec- 
tions pertain to general solutions of the stationary KdV hierarchy, we now concen- 
trate on the additional restriction that q be an elliptic function and hence return to 
our main subject, elliptic finite-band potentials q for L — d 2 /dx 2 + q(x), or, equiva- 
lently, elliptic solutions of the sta tionar y KdV hierarchy. The remarkable fin ite-g ap 
example of the Lame potentials ( 2.6.1 ) due to Hermite [113 and Halphcn ||l09| in 
the last century, brought back into the limelight by Ince [118 around 1940, remained 

1 3 

the only explicit elliptic finite-gap example until the KdV flow q t — —q X xx + -^QQx, 

with initial condition q(x,0) = — 6p(x + W3), was explicitly integrated by Dub rovin 
and Novikov §7| in 1975 (see also Enol'skii Its and Enol'skii JT||) and 

found to be of the type 

3 

(2.8.1) q(x,t) = -2^2p(x-Xj{t)) 

3=1 

for appropriate {xj(t)}i<j<3- Due to the unitary evolution operator U n (t) con- 
structed with the help of P2n+i(t) via U n .t{t) = P2n+i{t)U n (t) , all potentials q(x,t) 
in ( 2.8.1 ) are isospectral to q(x, 0) = — 6p(x + W3). 

In 1977, Airault, McKean and Moser, in their seminal paper ||, presented the 
first systematic study of the isospectral torus Iu(qo) of real- valued smooth potentials 
qa(x) of the type 

M 

(2.8.2) q (x) =-2J2p{x-Xj) 

i=i 

with a finite-gap spectrum. Among a variety of results they proved that any ele- 
ment q of Ir(<7o) is an elliptic function of the type ( 2.8.2 ) (with different Xj), with 
M constant throughout Iwi(qo) and dim/R(go) < M. In particular, if qo evolves 
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according to any equation of the KdV hierarchy it remains an elliptic finite-gap 
potential. The potential (2.8.2) is intimately connected with completely integrable 
many-body systems of the Calogero-Moser-type , ]166|] (s ee also Bennequin [H , 
Birnir |26| , Calogero |36) , Chudnovsky and Chudnovsky |39| , jy] , Olshanetsky and 



Perelomov [ 1 71 1 , and Ruijsenaars |189|). This connection with integrable particle 



systems was subsequently exploited by Krichever [ 136 1 in his construction of elliptic 
algebro-geometric solutions of the Kadomtsev-Petviashvili equation. In the KdV 
context of (2.8.2), Krichever's approach relies on the ansatz 



(2.8.3) 



3=1 



for the Floquet solutions of Lq = d 2 jdx 2 + qo{x), where 

a{x)cx{-p) 

(assuming for simplicity the generic case Xj ^ x^ (mod A) for j ^ k). Applying 
Lq to ( p. 8. 3 ) then yields an M-sheeted covering of the torus associated with the 
fundamental periods 2u>i , 2<x>3 and hence a description of t he u nderlying algebraic 
curve. (We will briefly comment on this ansatz in Remark 3.9.) 

The next breakthrough occurred in 1988 when Verdier [ 230 1 published new ex- 
plicit examples of elliptic finite-gap potentials. Verdicr's examples spurr ed a flurry 
of activities and inspired Belokolos and Enol'skii p9|, fl20f, Smirnov f|l96|, and 



subsequently Taimanov [208] and Kostov and Enol'skii [ L 30 1 to find further such 
examples by combining the reduction process of Abelian integrals to elliptic inte- 
grals (see Babich, Bobenko and Matveev |p|, jl4|| , Belokolos, Bobenko, Enol'skii, 
Its, and Matveev 0, Ch. 7, and Belokolos, Bobe nko, M ateev , and Enol'skii [p^| ) 
with the aforementioned techniques of Krichever [ 136 1 , [ 137 1 . This de velo p ment 
finally culminated in a series of recent results of Trcibich and Verdier |225|-[ p25| 
where it was shown that a general complex- valued potential of the form 



q{x) 



4 



(ll>2 = tjj\ + u>3, o->4 = 0) is a finite-gap potential if and only if dj /2 are triangular 
numbers, that is, if and only if 

dj = Sj(sj + 1) for some Sj € Z, 1 < j < 4. 

We shall from now on refer to potentials of the type 



q(x) 



i=i 



j(sj + l)p(x - ujj), Sj G Z, 1 < j < 4 



as Treibich- Verdier potentials. 
Colombo, Pirola, and Previato |4 



The methods of Treibich and Verdier (see also 
[185 1, Previato and Verdier 



Previato |184 



187 ) are based on hyperelliptic tangent covers of the torus C/A (A being the 



period lattice generated by 2wi and 2^3). 

The state of the art of elliptic finite-gap solutions up to 1993 was recently re- 
viewed in a special issue of Acta Applicandae Math., see, for instance, Belokolos and 
Enol'skii Q, Enol'skii and Kostov §§, Krichever [^39|, Smirnov (^98), Taimanov 
209 , and Treibich [ 22 1 1 . For more recent results see Eilbeck and Enol'skii J59[, 
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|60| , ]64]] and Smirnov [ 199 1 , [ 202 1 . Since a complete characterization of all elliptic 
finite-gap solutions of the stationary KdV hierarchy was still open at that time, 
we developed a new approach to this characterization problem to be described in 
Sections || and ||. As alluded to at the end of our introduction, Caloger-Moser-type 
models are again an intensive object of study. 

Since we will also discuss stationary rational solutions of the KdV hierarchy in 
Section U we should mention the case where p(x) in ( [2.8.2 ) degenerates into x~ 2 as 



discussed, for instance by Airault, McKean, and Moser |5j], Krichever ]135[ , [ 138 1 , 



Moser [166], [167|, Pelinovsky [179], and Shiota E9 



2.9. Linear Differential Equations in the Complex Domain. While all the 
developments described in previous subsections were in place around 1993, one final 
point, the connection of this subject to the classical area of differential equations 
in the complex domain, was made only around 1994 when we started to work on 
|98| . In the following we will remind the reader about this fundamental, but thus 
far missing piece, which plays a decisive role in the remainder of this review. 

In the late 1830's, Lame studied Laplace's equation in confocal coordinates. After 



some appropriate changes of variables this led to the differential equation (2.6.1), 
it's Weierstrass form. At the end of the last century, Lame's equation (especially in 
Jacobi's form) was studied intensively by Hermite [112], who obtained the general 
solution of the equation for integer values of n and any value of z. Picard proceeded 
to consider first general second-order equations fL80[, then n-th order equations 
see also Floquet 



181 



109 



J3|-|75J, Mittag-Leffler |63|, and Halphen g07|], |08f, 
), and finally first-order systems [182] whose coefficients are elliptic functions. 



Consider the differential equation 



(2.9.1) 



y '(x) = Q(x)y(x), 



where y(x) is C™-valued and Q is a n x n matrix whose entries are elliptic functions 
with a common period lattice spanned by the fundamental periods 2lu\ and 2^3. Let 
y be the space of solutions of (2.9.1) and denote the restriction of the translation 



operator y 1— > y(- + 2ujj) to y by Tj. Using this notation, Picard's theorem reads 

Ch. 15, Forsyth f|, 
Sect. 15.6, Krause 



as follows (se e also Akhiezer H, Sects. 58, 59, Burkardt [|3| 



Ch. IX, Gray (103 
132) , Vol. 2, Ch. 3 



Sect. 6.1, Halphen (109|, Ch. XIII, Ince (119 
and Picard p3|, Sect. III.V) 



Theorem 2.7. Assume that the first- order system (2.9.1) has a meromorphic fun- 
damental system of solutions. Then there exists at least one solution ip which 
is elliptic of the second kind, that is, the components of ip are meromorphic and 
Tjtjj = pjtjj for j — 1,3 for suitable constants pi, P3 €E C\{0}. // in addition, one of 
the operators T\ and T3 has distinct eigenvalues, then there exists a fundamental 
system of solutions of (2.9.1) which are elliptic of the second kind. 



The explicit Floquet-type structure of solutions of (2.9.1) in terms of a dou- 
bly periodic vector, powers of x, powers of the Weierstrass zeta-function, and an 
exponential contribution, has recently been determined in fl9l| . 

About the time Floquet, Fuchs, Hermite, Mittag-Leffler, and Picard (cf. the 
historical discussions in Gray |103| , Ch. VI) developed the theory of differential 
equations with elliptic coefficients, Floquet ]7l| , (72) also established his celebrated 
results for linear, homogeneous differential equations with simply-periodic coeffi- 
cients. From a historical perspective it is perhaps interesting to note that Floquet 
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assumed the coefficients, as well as the general solution of the equation, to be mero- 
morphic in order to arrive at the existence of periodic solutions of the second kind, 
that is, he obtained the precise analog of Picard's result. Only later was it realized 
that his theorem extends to continuous periodic coefficients on K without any refer- 
ence to meromorphic fundamental systems. The solutions Floquet called "periodic 
of the second kind", are today generally called Floquet solutions. 

Next, we mention another theorem with a similar flavor that concerns differen- 
tial equations with rational coefficients and meromorphic fundamental systems of 
solutions and hence is applicable to the study of rational algebro-geometric KdV 



potentials. In 1885, Halphen 1 108 published the following result 



Assume that the differential equation y( n >(x) + qi(x)y^ n ~ 1 ' (x) + 
= has rational coefficients bounded at infinity and a meromor- 
phic fundamental system of solutions. Then the general solution is of the form 
Ri(x) exp(Aix) + • • • + R n {x) exp(A„a;), where R\, . . . , R n are rational functions. 



Theorem 2.8 

h q n (x)y(x 



It should be emphasized that the principal hypothesis in Theorems [2. 7] and 2 
the existence of a meromorphic fundamental system of solutions, can be verified 
in a straightforward manner by applying the Frobenius method (see, e.g., Cod- 
dington and Levinson @, Sect. 4.8, Forsyth |7(|, Ch. Ill, Hille |ju||, Ch. 9, Incc 
|TT9| , Ch. XVI, and Whittaker and Watson |235fl , Ch. X) to each pole of q in the 
fundamental period parallelogram. 

Finally we mention an observation made by Appell . Let y\ (x) and yi (x) be 
linearly independent solutions of y"(x) + u{x)y(x) — 0. Then y\{x) 2 , yi(x)y2{x), 
and y2(x) 2 are linearly independent solutions of 



(2.9.2) 



w (x) + 4m(x)w (x) + 2u (x)w(x) = 0. 



This equation is easily integrated and yields 

(2.9.3) g'{x) 2 - 2g(x)g"(x) - Au{x)g{xf = W( yi ,y 2 ) 2 , 

where g(x) is the product of any two solutions y\{x) and y2(x) of y"(x)+u(x)y(x) = 
and W(yi,y2) denotes their (x- independent) Wronskian. This innocent looking 
fact will be of great importance in our analysis later. In fact, a comparison of 
equations ( 2.3.4 ) and ( 2.9.2 ) reveals another connection between the KdV hierarchy 
and the Schrodinger-type equation y" (x)+u(x)y(x) = 0. Moreover, since the formal 
Green's function G(z,x, ,x') of d 2 /dx 2 + q(x) on the diagonal x = x' is of the 



type yi(z, x)y2(z, x)/W(yi(z), y2(z)), ( 2.9.3[ ), with u(x) = q(x) — z, is equivalent to 
the well-known universal nonlinear second-order differential equation satisfied by 
G(z,x,x) (see, also Gelfand and Dickey, |^3|, f85|). 

It should be noted that Drach (49|, Q, |51j (see also used ( |2 . 9 . 2| ) to derive a 
class of completely integrable systems now known as the stationary KdV hierarchy 
as early as 1918/19. It appears he was the first to make the explicit connection 
between completely integrable systems and spectral theory. More than 55 years 
later, Gelfand and Dickey l83[, p5[ also based some of their celebrated work on the 



KdV hierarchy on (2.9.2) 



3. Algebro-Geometric, and Especially, Elliptic KdV Potentials. 

Definition 3.1. Suppose q is meromorphic and let L be the differential expres- 
sion L = d 2 /dx 2 + q(x). Then q is called an algebro-geometric KdV potential (or 
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simply algebro-geometric) if q is a solution of some equation of the stationary KdV 
hierarchy. 

Equivalently, we could define the meromorphic function q as algebro-geometric 
if any one of the following three conditions is satisfied. 

(1) There exists an odd-order differential expression P such that [P, L] = (ac- 
cording to the result of Burchnall and Chaundy) . 

(2) There exists an ordinary differential expression P of odd order and a polyno- 
mial R such that P 2 = R{L). 

(3) There exists a function F : C 2 — > Coo, which is a polynomial in the first 
variable, meromorphic in the second, and which satisfies F'" '(z, x) + 4(q(x) — 



z)F'(z,x) +2q'(x)F(z,x) = (cf. equation Q2.3.4Q ). 

It can be shown (see Theorem 6.10 by Segal and Wilson Jl94[ ) that any solution 
of any of the stationary KdV equations is necessarily meromorphic. Hence the 



assumption that q is meromorphic is actually redundant in Definition 3.1 



3.1. Periodic KdV Potentials. If q is real-valued, locally integrable, and pe- 
riodic we obtain the equivalence of the following statements from the works of 
Birkhoff and Hamel described in Section |2.6| . 

(i) q is algebro-geometric with P 2 = JJ^oi^ ~ ^j')- 

(ii) q is finite-band with spectrum the union of (— oo, A2J and n compact bands 
[Aaj+i, A 2 j], 3 = 0,...,n- 1. 

(iii) Ly — zy has two linearly independent Floquet solutions for all z £ C with 
the exception of the 2n + 1 values Ao, A2 n - 

(iv) q has n movable Dirichlet eigenvalues, precisely one in each of the closures of 
the spectral gaps (A2j, A2j-i), j = 1, n. 

As mentioned previously, the classical example for finite-band potentials are the 
Lame potentials q(x) = —s(s + l)p(x + c;wi,u>s), n £ N, with c — u>3 purely 
imaginary and u>\ real. However, Lame potentials are algebro-geometric for general 
choices of the half-periods as well as for general choices of c £ C. This suggests the 
study of complex-valued potentials with inverse square singularities as in and 



233 , and we will subsequently report on some of these results. 

Let q be a complex- valued, periodic function with period Q > 0, which is locally 
integrable on R\E, where E G K is a discrete set (i.e., a set without finite accumu- 
lation points). Moreover, q is assumed to be meromorphic near each £ G E with 
principal part — s(s+l)/(x— £) 2 , where s — s(£) £ N. Then one can define unique so- 
lutions of initial value problems of the differential equation y" (x)+q(x)y(x) = zy(x) 
on R\E (with initial conditions at xq £ R\E) by analytic continuation around the 
singularities. Even though the potential is no longer continuous, Floquet theory 
(see Section remains essentially unchanged. 

To apply Floquet theory we first introduce a basis in y(z), the space of solutions 
of y"{x) + q{x)y{x) — zy(x). Let c(z, -,xq), s(z, -,xq) £ y(z) be defined by the 
initial conditions c(z, xq, xq) = s'(z, xq, xq) — 1 and c'(z,xq,xq) = s(z, xq,xq) = 
(prime denoting the derivative with respect to the second variable). Let p± denote 
the Floquet multipliers of y"{x) + q{x)y(x) = zy(x), that is, the eigenvalues of 
the translation operator 7q(z). An important role is played by tr(7n(z)), which is 
sometimes called the Floquet discriminant and which, in our basis, is given by 

tr(7fl(z)) = C(Z, Xq + Q, Xq) + s'(z, Xq + f2, Xq). 
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Since the trace is of course independent of the chosen basis in y(z), the dependence 
of the right-hand side on xo is only apparent. The Floquet solutions, may be 
expressed as 

f±{z,x ,x) = s(z,x + Q.,x Q )c(z,x Q ,x) + {p± - c(z,x + Q, x ))s(z, x , x), 
and their Wronskian is given by 

W(f + (z,;X ),f-(z,;Xo)) = -s(z, x a + Q, x ) V(tr(T n (z))) 2 - 4. 
Next, consider the function 

/„ -, t \ f+{z,x,xo)f-{z,x,x ) 

3-1.1 g{z,X) = r - F - ? rr. 

W{f + {z,-,xo),f-{z,-,x )) 

As our notation for g(z,x) suggests, the dependence of this function on xq is only 
app arent since f±(z, •, x\) are just multiples of f±(z, •, xo) and the right-hand side 
of (|3.1.lD is independent of normalization. In particular, we may replace xq by x 



in (3.1.1) to obtain 

-s(z, x + 0, x) 



g{z,x) 



v /(tr(T n (z))) 2 -4 

The function s(-,x + f2, x) is an entire function of order of growth 1/2. The zeros 
of s(-,x + fl, x) are the Dirichlet eigenvalues of d 2 jdx 2 + q on the interval [x, x + Q] 
and their order, which we denote by d(z, x), is the algebraic multiplicity of the 
corresponding Dirichlet eigenvalue (we set d(z,x) = if s(z,x + Q,x) ^ 0). We 
also introduce ck(z) = mxa{d(z, x) \ x € R\S} and d m {z,x) = d{z,x) — di{z), the 
immovable and movable parts of d{z, x), respectively. The quantity J2zec x )i 
which is independent of x, is called the number of movable Dirichlet eigenvalues. 
Using Hadamard's factorization theorem we write g(z,x) — F(z, x)D(z) collecting 
in F the factors depending on x and in D the factors independent of x. Then the 
multiplicity of a zero z of D is just di(z), while the multiplicity of a zero z of F(-, x) 



is d m {z,x). One then obtains from Appell's equation (2.9.3) that 
(3.1.2) 

F'(z, xf - 2F(z, x)F"(z, x) - A(q - z)F(z, xf = ((tr(T^(z))) 2 - A)/D{zf. 

Recall that a zero z of (tr(Tn(z))) 2 — 4 is an (anti-)periodic eigenvalue whose mul- 
tiplicity we denote by p(z). Since the left-hand side of equation ( p. 1.2 ) is entire as 
a function of z we obtain the following result. 

Theorem 3.2. There exists an entire function R such that (tr(7n(z))) 2 — 4 = R{z) 
D(z) 2 . In particular, p(z) — 2di(z) > for every zeC. 

There are, at most, countably many points where p(z) > since these points are 
isolated. Therefore, there are at most countably many points where p(z) — 2di(z) > 
0. These include all algebraically simple (anti-)periodic eigenvalues (where p(z) = 1) 
but may well include other points too. 

We call a Floquet solution ip(zo, •) of y"{x) + q{x)y(x) = z y(x) regular if there 
exist Floquet solutions ip(z, •) of y" (x)+q(x)y( x) — zy(x), which converge pointwise 
to ip( z 0i ') as z tends to zq. It was shown in ]233| that the set of regular Floquet 
solutions forms a line bundle on the topological space Mp obtained from the curve 
p 2 — (tr(Tn(z)))p+ 1 = by desingularization at all points where p{z) = 2d%{z) > 0. 
The space Mf can be viewed as a double cover of the complex plane branched 
precisely at all points z where p(z) — 2di{z) > 0. In particular, the equation y"(x) + 
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q(x)y(x) = zy(x) has two linearly independent regular Floquet solutions if and 
only if p(z) — 2di(z) = 0. In other words, p(z) — 2di(z) > indicates a defect in the 
structure of regular Floquet solutions. The number def(g) = ^2 ztE c(p{z) — 2di(z)) — 
deg(R), which is a positive integer or infinity, will therefore subsequently be called 
the Floquet defect of q. When q is real- valued and nonsingular, then p(z) — 2di{z) 
is always zero, except when p(z) — 1, which forces dt(z) = 0. In this case dcf(<7) 
counts the number of points z G C where only one linearly independent Floquet 
solution exists. In general, however, y" (x) +q(x)y(x) — zy(x) may have two linearly 
independent Floquet solutions with only one being regular. 

Theorem 3.3. The following statements are equivalent: 

(1) def(q) =2n+l. 

(2) There are n movable Dirichlet eigenvalues. 

(3) There exists a differential expression P271+1 of order 2n+l but none of smaller 
odd order commuting with L. This differential expression satisfies -Rfn+l = 
i?2n+i(z) = Ilzec^ — z) p ^~ 2di ^ and hence q is algebro- geometric. 

Sketch of proof. Generally p(z) < 2 and d(z, xq) < 1 for all suitably large z. More- 
over, asymptotically, any Dirichlet eigenvalue is close to two (anti-)periodic eigen- 
values. Therefore, def(g) = 2n + 1 implies d m {z, xo) is different from zero for only 
finitely many z G C and J2 z ec dm(z, x), which equals the degree of F(-,xq), must 



be finite. Equation (3.1.2) then yields degF(-,a;o) = n. Hence (1) implies (2). The 
converse of this follows immedi ately from equation ( 3.1.2| ). 



Differentiation of equation ( 3.1.2|) shows that the third criterion after Defini- 



tion is satisfied when F(-,xq) is a polynomial. Hence (2) implies (3). We 
remark here that i?2n+i is a constant multiple of R. To prove that (3) implies (2) 
one has to show that the zeros of the function F n defined by P2n+i are precisely 
the mova ble Dirichlet eigenvalues of L. This follows from applying Pm+\ given 
by ( 2.3.2] ), successively to the generalized Dirichlet eigenfunctions (i.e., the eigen- 



functions corresponding to the algebraic eigenspace) associated with the movable 
Dirichlet eigenvalues of L (cf. |233| for more details). □ 

Remark 3.4. The quantity p(z) — 2di(z) proves to be of utmost importance. Deter- 
mining the points where it takes on positive values, this quantity then governs the 
structure of the line bundle of regular Floquet solutions and determines the entire 
function R(z) = ((tr(7n(z))) 2 — 4)/D(z) 2 , which defines the algebraic curve associ- 
ated with q if q is algebro-geometric and R is a polynomial. When q is real-valued 
and locally integrable on R, then geometric and algebraic multiplicities of Dirichlet 
and (anti-)periodic boundary value problems coincide. Therefore, < p(z) < 2 
and < d(z,xo) < 1. Thus, di(z) — 1 enforces p(z) = 2. In addition, if p(z) = 2, 
then all solutions are (anti-)periodic, implying di(z) = 1. In this case the questions, 
"When is p{z) — 2di(z) > 0?" and "When is p(z) — 1?", are equivalent. Hence, in 
determining the edges of the spectral gaps, the role played by the Dirichlet eigen- 
values and, in particular, the distinction between movable and immovable Dirichlet 
eigenvalues, is secondary in the case of real-valued locally integrable potentials. 

Remark 3.5. One may also show that every algebro-geometric potential is a finite- 
band potential, that is, the conditional stability set (which coincides with the spec- 
trum of H when no singularities are present) consists of finitely many regular ana- 
lytic arcs. 
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When q is real and has no singularities, the converse is also true (Dubrovin f33|]). 
However, in general this is not the case as the following example shows. Let q(x) = 
e 2lx . Then a fundamental system of solutions of y"(x)+q(x)y(x) = — v 2 y(x) is given 
in terms of Bessel functions (cf. j|, Ch. 9) by yi(x) — J v {ie lx ), y2{x) — Y v (i% %x ). 
Note that y\ is always a Floquet solution with multiplier e vv% . Hence z = — v 2 is 
in the conditional stability set if and only if v € R. Consequently, S(L) = (—00, 0]. 
However, s(i> 2 ,xq + n 7 x n ) = ■KJ u (ie tXo )J- v (ie ix °), which is entire as a function of 
xq. Hence di(—v 2 ) = for all v € C, that is, every Dirichlet eigenvalue is movable. 
Thus, def(g) = 00, and hence q is not algebro-geometric. More general examples 
of this type have been studied systematically by Gasymov p0| , j8l| , Guillemin and 



Uribe |10§, and Pastur and Tkachenko pT^ , p75 |. 
3.2. Picard-KdV Potentials. 

Definition 3.6. Let q be an elliptic function. Then q is called a Picard-KdV 
■potential (or simply a Picard potential) if the equation y"(x) + q(x)y(x) = zy{x) 
has a meromorphic fundamental system of solutions with respect to x for all values 
of the spectral parameter z G C. 

Theorem 3.7. If q is a Picard potential then it may be represented as 

m 

(3.2.1) q(x) =C-J2 + 1)p0& - *>i) 

3=1 

for suitable integers si,...,s m and complex numbers C, bi,...,b m , where the bj are 
pairwise distinct mod A. 

Sketch of proof. Every singularity of y"(x) + q(x)y(x) = zy(x) must be a regular 
singular point with integer indices. From the partial fraction expansion for elliptic 



functions (Theorem 2.2) one obtains 



q(x) =C- 5>,(s, + l)p(x - bj) + BjCix - bj)). 

3=1 

In a vicinity of bj there is a solution of the form 

00 

1>(x) - (x-bj)- s ^p k (x-bj) k , 

k=0 

where (3q = 1. Next we use the Frobenius method to show that Bj = 0. Let 

v ■>' ■> fc=0 

and insert tp into the differential equation y" + (q — z)y — to get 

=f(-Sj)(x - bj)- s ' + {/(l - Sj)f3! + G 1 }(x - b,) 1 -^ 

+ ■■■ + {/(* - Sj)p k + G k }(z - bj) k - s * +■■■ , 

where f(l) = (£+ Sj)(t - sj - 1) and 

G k = BjP k ^ + (Cj, - z)/3 k -2 + Gj,lP k - 3 + ■■■ + Cj,k-20Q- 

Now f(—Sj) = and j3\, (32 Sj may be determined recursively so that the coeffi- 
cients of (x — bj) 1 ^ Sj , (x — bj) Sj vanish. But since f(sj + 1) = 0, the coefficient 

b 3 r : 



of (x — bj) Sj+1 is just G2SA+1, which therefore must vanish for all z £ C. On the 
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other hand, if Bj ^ 0, one can show by induction that G2 Sj +i is a polynomial in z 
of degree Sj. This contradiction completes the proof. □ 

If q is a Picard potential then, by Picard's theorem (Theorem ^?]), the equation 
y"(x) + q(x)y(x) = zy(x) has at least one solution which is elliptic of the second 
kind. Using Theorem [^4| and the special structure of the Picard potential (3.2.1) 
this solution may be represented as 

i~n=i a ( x ~ a j( z )) 

(3-2.2) = -FFm 7 TT7- ex P(^a( z )2:), 

where s ~ - —1 Sj and a(z) — (ai(z), . . . , a s (z)). At bj the function ip a has a pole 
of order Sj or a zero of order + 1. For later notational purposes we allow for 
Sj = 0, in which case tp a has either no pole and no zero or a simple zero. For 
subsequent use we define 

Mi = {j e {1, ...,m} | ordf^a) = sj}, 

M 2 = {j G {1, m} | ordfe. (V-a) - Sj + 1}. 

The function tp a is a solution of y" {x) + q(x)y(x) = zy(x) if and only if 

s m 

(3.2.3) K + ^ " °i) - £ - 6 i) = °> 

3=1 3=1 

s m 

(3.2.4) z = C - ^(1 - 2s r )p(b r - aj) - ^ Sj{ Sj + 2s r )p(b r ~ bj), 

3=1 3=1 

where r is chosen such that s r ^ and ordb r (ipa) — —s r which is always possible. 
In the case of a Lame potential these conditions are recorded, for instance, by 
Burkhardt |33j. The subscript a in A a expresses the dependence of A a on a = 
(ai, a s ), which in turn depends on z. 

In |)6| we have developed a method for even Picard potentials (i.e., potentials q 
satisfying q(xo + x)) — <l{ x o ~ x) for some :ro G C) to determine all points where 
two regular Floquet solutions fail to exist. (For simplicity we will assume xq = 
from now on.) First note that an even Picard potential is of the form 

(3.2.5) 

4 m 

q(x) = C - ^s fc (s fc + l)p(x - uj k ) - ^r fe (r fe + l)[p(x - b k ) + p(x + b k )}, 

k=l k=l 

where the bk are pairwise distinct and different from half-periods, the Sk are non- 
negative and the are positive. If m = 0, that is, if 

4 

q(x) = C - ^ s k{$k + l)p(z - Wfc), 

fc=l 

then q is called a Treibich-Verdier potential following the work of Verdier J230| 
and Treibich and Verdier |21S|— [227|. If in addition, only one of the numbers s k is 
different from zero, then q is a Lame potential. In order to make use of previous 
results we will adopt the following notation: m = 2fh + 4, bk+fh = —bk, Sfc+m = r k 
for k = 1, ...,m, and bk+2m = <^k, Sfc+2m = s k for k = 1, ...,4. Let s be the number 
of Oj's which do not appear in {b\, b m }. 
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If q is an even Picard potential of the type (3.2.5), and iiipa given by ( 3.2.2| ) is a 
solution of the differential equation y" (x)+q(x)y(x) = zy{x), then so is the function 
ip-a, which is obtained by replacing every aj with —aj in ( |3 . 2 . 2| ) and ( |3.2.3| ), since 
ijj_ a {x) = (-l) sl + S2+s '^ a (-x). 

Next we compute the Wronskian of the two solutions ip a and tp—a- One obtains 
an expression which involves x, but since the Wronskian does not depend on x, one 
may evaluate it anywhere. Choosing x — ag for any ag which does not appear in 
...,b m }, one finds that the Wronskian is a nonzero multiple of 



a{2ai) 



a{ag - b r )a(ag + b r ) 

Since ag is different from all the bj and, in particular, different from the half-periods, 
one infers a(2ag) ^ 0. Also ag ^ if k ^ I. Therefore we find that the Wronskian 
is zero if and only if u(ag + aj) = for some j S {1, ...,£ — 1,£+ 1, s} and hence 
aj = -a? (modA). In particular, we find that the number s is even and define 
d = s/2. 

Choosing now x — bg for any £ £E M2 the Wronskian can be written as a nonzero 
multiple of 



<y{2bg) 



2se+ i U 3 eM 2 o(bg - bjf^a{bg + b,j)^+ l 



a(bg-b r )a(bg + b r ) I\7=i a ( b £ ~ ^) 2Sj 

which is zero if and only if bg is a half-period or if there is a j € Mi such that 
bj = —bg (mod A). 

In summary we have found the following: if ip a and tp- a are linearly dependent 
solutions of tj)" +qip = zip, then some of the numbers ai, a s may be half-periods 
while all others appear in pairs (aj, ag,.) with ag j = —aj. Moreover, if aj is equal 
to a half-period u>k, which is a pole of q of the form — rfe(r^ + l)/(z — Wfe) 2 , then 
there are exactly 2rk + 1 of the ag which are equal to this half-period. If aj equals 
a pole bg of the form —sg(sg + \)/(z— bg) 2 , where bg is not a half-period, then there 
are exactly 2sg + 1 of the a m which are equal to this pole and exactly 2sg + 1 other 
a m 's which are equal to the pole — bg. 

This information is now being used to rewrite the solution -ipa of ip" + qip = zip for 
those values of the spectral parameter z where W(ip a , ip-a) — as a product of two 
functions. The first function is fixed, depending only on the poles of the potential q, 
on the half-periods, and the exponents associated with these. The second function 
is a polynomial in p(x), whose coefficients depend on those aj which are neither 
half-periods nor poles of q and which are as yet undetermined. According to the 
above argument there must be an even number, 2d, of those, and half of them are 
just negatives of the other half. 

Hence we define tg — OTdb e (ip a ) for £ = 1, ...,2m -I- 4 and obtain t^ a {x) — 
f(x)Q(p(x)) 7 where 

f(x) = e A °* ( H a(x - 0*)**+** ] ( n a(x - bg)A a(x) 



\k=l 

d 



2,1 



Q(p(x)) = Y[(p(x) - p(flj)) = ^c,jp(xy 
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Here we used the fact that a(z — a,j)cr(z + a,j) — —<r(z) 2 a(aj) 2 (p(z) — p(a,j)). 
Moreover, we dropped the non-zero constant factor (— l) d IIjLi a ( a j) 2 - This yields 

for suitable constants Sj t k and Pj.k depending on q and the numbers tg. Let S and 
T be the matrices with entries Sj,h and Tj t k- 

Hence we obtain a solution ip a of the equation y" (x)+q(x)y(x) — zy{x) satisfying 
W(ip a , i>—a) — if and only if 

S"/ = 27 and T7 = 0, 

where 7 = (c , ■ ■■ 1 c d ) T ■ 

For any given even Picard potential there are several (but finitely many) choices 
to distribute some (or all) of the parameters a±, a s among the half-periods and/or 
poles of q. Accordingly, there are several (but finitely many) of the above described 
constraint eigenvalue problems to solve in order to find all the values of the spectral 
parameter z € <C where W(ip a ,ip— a ) = 0. In each case there are only finitely many 
eigenvalues of the associated matrix S, some (or possibly all) of which may be in 
contradiction to the constraint T7 = 0. Thus, we proved the following result. 

Theorem 3.8. Let q be an even Picard potential. Then, for every complex number 
z, the differential equation y"(x) + q{x)y{x) = zy{x) has a solution tp a of the form 
(HJ), that is, a solution which is elliptic of the second kind. Similarly, the function 
ip- a is a solution of the same equation (for the same value of z) and also elliptic 
of the second kind. With respect to any period of q, the functions ip a and ip- a are 
regular Floquet solutions. For all but a finite number of values of z £ C these two 
solutions are linearly independent and therefore p(z) — 2di(z) — 0. Hence def(g) < 
00, that is, q is algebro-geometric. 

In the case of a Treibich-Verdier potential 

4 

(3.2.6) q(x)=-Y,s j (s j + l)p{x-u> j ), Sj G N U {0}, j = 1, . . . ,4, 

the matrix T is absent and therefore one only has to find the eigenvalues of S 
for all possible choices of the numbers fy. This was performed in |Mj for Lame 
potentials and in (95| for Treibich-Verdier potentials. The arithmetic genus of 
the curve P 2 = R{L) associated with q is given in the following table, where 
s = si + S2 + S3 + S4 (and, without loss of generality, si > S2 > S3 > S4 > 0). 



s 




f/= of finite branch points 


genus 


even 


S2 + S3 < Si + S 4 


2si + 1 


Sl 


even 


S2 + S3 > Si + Si 


Sl + s 2 + S3 — S4 + 1 


f -84 


odd 


S2 + S3 + S 4 < Si 


2si + 1 


Sl 


odd 


S2 + S3 + S 4 > Si 


si + s 2 + s 3 + s 4 + 2 


s+4 

2 



Shortly after this table appeared in our paper [ )95[ , Armando Treibich kindly sent 
us the following formula for the arithmetic genus of the Treibich-Verdier curve. It 
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is given by 

max{2 Sl , s + 1 - (1 + (-l) s )(s 4 + (l/2))}/2 
in agreement with our column on the right. 



Remark 3.9. We briefly return to the Picard potential q(x) in ( 3.2.1 ). Assuming 
first the special case where Sj = 1, j = 1, . . . , m, that is, 

m 

q(x) = P(« ~ bj), b k ^ bt (mod A) for k f I , 

3=1 

Theorem 2J5 yields Krichever's ansatz ( 2.8. 3| ) in [ 136 1 , that is, 

a(x — cij(z)) 



(3.2.7) ip a ( z )(x) = exp(A a(z) .T) 
where 



_ e n{z)x 



$(x, p) = e c(p)a; CT(a; - p)a{x)- 1 a{- p)~ 



and where k, Ay and pj are suitably chosen. 
In the general case, where 



(3.2.8) 



q(x) 



Sj(sj + l)p(x - bj), b k be (mod A) for k ^ £, 

3 = 1 



one can use an extension of Theorem 2.5 (using, e.g., l'Hospital's rule as several 
bj's are confluenting into a bj and converting 6-derivatives into an x-derivative) to 
arrive at the corresponding analog 



(3.2.9) 



i>a{z){x) = exp(A a(z) a;) 



Yl7=i a ( x - b i) Sj 



j=l k=0 



3 = 1 



of ( 3.2.7 ). The extended ansatz ( 3.2.9 ) was recently used by Enol'skii and Eilbeck 
59 1 . |64[, Enol'skii and Kostov |65|, |13C|, in the context of Treibi ch-V erdier poten- 
tials. The special case m = 1, Si = 2 can be found in H crmitc [113, p. 374-377 
(see also For syth |7q| , p. 475-476), the general case ( 3.2.9 ) is discussed in Krause's 
monograph (l32), Vol. 1, p. 292-296, Vol. 2, p. 183, 259-264. 

We emphasize that Picard's theorem, Theorem 2.7, yields the er-function rep- 
resentations of ipg( z) (x) in ( 3.2.7 ) and ( |3.2.£ ) and hence Krichever's ansatz ( 2.8.3 ) 
and its extension ( 3.2.9| ). In particular, an alternative characterization ofah sta- 
tionary elliptic KdV and AKNS potentials to the one provided in Theorem |3.12 and 
Theorem 4.8, respectively, can be based on these observations (cf. 1 101 ] ) . In the 
remainder of this review we do not pursue this avenue but stress Floquet-theoretic 
and Green's function methods instead. 

3.3. Necessary Conditions for a KdV Potential to be Algebro-Geometric. 

Theorem 3.10. Suppose q is an algebro- geometric KdV potential. Then (Ly)(x) = 
y"(x) + q(x)y(x) = zy(x) has a meromorphic fundamental system of solutions with 
respect to x for all values of the spectral parameter z G C. 
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Sketch of proof. First one shows that every pole of q is of second order. This follows 
since only in this case there is a balance between the growth of the order of the 



poles of /•" and 4<?/j + 2q' fj considering the recursion relation ( |2.3.3| ). Such a 
balance is necessary since eventually all fj vanish. Next one shows by a similar 
argument that the leading term of the Laurent expansion of q about a pole Xq 
must be — s(s + l)/(x — xq) 2 for some integer s. This guarantees that the exponents 
of the singularity xq of the differential equation y"(x) + q(x)y{x) = zy(x) are 
integers. Finally one has to show the absence of logarithmic terms in the solutions 
of y"(x)+q(x)y(x) = zy(x) which follows from a careful analysis using the Frobenius 



method. More details are provided in [234]. □ 



We note that r-function results in Segal and Wilson 1 194 imply Theorem |3 . 10 
(In the case of nonsingular curves K. n this simply follows from the standard thcta- 
function representation of the Baker- Akhiezer function.) While Segal and Wilson 
rely on loop group techniques (and study the Gelfand-Dickey hierarchy), the above 
proof represents a completely elementary alternative for the KdV hierarchy. 

3.4. Sufficient Conditions for a KdV Potential to be Algebro-Geometric 
Characterizations of Elliptic, Simply Periodic, and Rational Stationary 
KdV Solutions. 

Theorem 3.11. q is an algebro- geometric potential if, for each z G C, y"(x) + 
q{x)y{x) = zy{x) has a meromorphic fundamental system of solutions with respect 
to x and if one of the following three conditions is satisfied. 

(i) q is rational and bounded near infinity. 

(ii) q is simply periodic with period VL and there exists a positive number R such 
that q is bounded in {x S C | | Im(a;/51)| > R}. 

(hi) q is elliptic. 

Sketch of proof. We will first sketch the proof in |98| which treats the case where q 
is elliptic and prominently uses the double periodicity of q. 

Let q be an elliptic function with fundamental periods 2ui and 2W3. Assume, 
without loss of generality, that Im^/wi) > 0. Introduce tj — LUj/\u>j\ and define 

qj(x) = t%q(tjx). 

The transformation £ — tjX, ip(t;) = w(x) transforms ip" + qip = zip into w" + qjW = 
t 2 zw. 

Note that qj has the period 2\ojj\. With the aid of Rouche's theorem, one may de- 
termine the asymptotic distribution of the (anti-)periodic eigenvalues of d 2 jdx 2 + qj 
and verify that they all stay close to the real axis while their re al pa rts tend to — oo. 



(Never mind that qj may have inverse square singularities, see [232].) Equivalently, 
all 2tjj-(anti-)periodic eigenvalues of d 2 jdx 2 + q lie in the half-strip £j given by 

Ej = {z e C 1 1 Im(^ 2 z )| < C, Rc(t 2 z) < C*}, j = 1, 3 

for some constant C > 0. Note that the angle between the strips Si and S3 is 
positive and less than 2ir and therefore, they intersect only in a finite part of 
the complex plane. Hence, for every point z outside a sufficiently large compact 
disk, at least one of the translation operators T\ or T% has distinct eigenvalues. 
Assuming, according to our hypothesis, that all solutions of y"{x) + q(x)y(x) = 
zyix) are meromorphic, we may invoke Picard's theorem to conclude that there 
are two linearly independent solutions which are elliptic of the second kind for 
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any z outside that disk. These solutions are Floquet solutions with respect to 
any period of q and since the disk is compact, we have shown that there are at 
most finitely many points which lack two linearly independent Floquet solutions. 



Applying Theorem 3.3 then proves that q is algebro-geometric. 

Next we sketch the proof for the case when q is rational. That proof can easily be 
extended for the simply periodic and elliptic potentials. In the latter case one uses 
the algebraic properties of elliptic functions rather than their double periodicity. 



For more details see [234 



Suppose that q is rational and bounded at infinity. Let z = lim^—Kx, q{x). From 
Halphen's theorem (Theorem |2.§| ) we obtain for z ^ z Q the existence of linearly 
independent solutions 



y±(z,x) = R±(z, x) exp(±V z — z^x), 

where R±(z, •) are rational functions. The poles of R±{z, •) are determined as the 
singular points of the differential equation and their orders as the exponents of the 
corresponding singularities. Next define the function g(z,x) = y + {z 1 x)y-(z 1 x). 
Then there exists a polynomial v in x such that 

d 

v{x) 2 g(z,x) = ^2c j (z)x j . 

Next note that the functions v 2 g(z,x), v 3 g'(z,x), v 4 g"(z,x), and v 5 g'"(z,x) are 
polynomials in x whose coefficients are homogeneous polynomials of degree one in 
Co, Cd- Also v 2 q and v 3 q' are polynomials. Hence v 5 (g'" + 4(g — z)g' + 2q'g) is 
also a polynomial in x, whose coefficients are homogeneous polynomials of degree 
one in cq, . . . , cj. The coefficients of the polynomials q in this last expression are 
polynomials in z of degree at most one, that is, 

N d 

(3.4.1) v 5 (g"' + A(q - z)g' + 2q'g) = £X>i,/ + Pjrfctx* 

j=0 £=0 

for suitable numbers N, (Xj,e, and (3j^ which only depend on q. From Appell's 



equation (2.9.2) it follows upon differentiation that the expression (3.4.1) vanishes 
identically. This gives rise to a homogeneous system of N + 1 linear equations which 
has a nontrivial solution. Solving the system shows now that the coefficients q are 
rational functions of z. Therefore 

F(z,x) 

g{z,x) = — — , 

7 (z) 

where F{-, x) and 7 are polynomials and F(z, ■) is a rational function. Hence q is 



algebro-geometric by item (3) following Definition 3.1. □ 



Combining the results of this subsection and its preceding one yields the following 
explicit characterization of all elliptic algebro-geometric KdV potentials (a problem 



posed, for instance, by Novikov, Mana kov, Pitaevskii, and Zakharov [170|, p. 152) 



originally proven in |98| (see also |97[ , |100| ) ■ 



Theorem 3.12. Let q be an elliptic function. Then q is an elliptic algebro-geo- 
metric KdV potential if and only if it is a Picard potential. 

Similarly, these results characterize the stationary rational KdV potentials van- 
ishing at infinity studied, for instance, by Adler and Moser H, Airault, McKean, 
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and Moser M , C aloge ro |}5| , Q . Ch udnovsky and Chud novsky |3l| , [|d| , Grinevich 
|l05|, Krichever pi, Matveev fL5j|, and Sokolov p07|. 



4. Algebro-Geometric and Especially, Elliptic AKNS Potentials. 

Most of the results in this section are taken from |pjj|| . Since we take here a 
slightly different point of departure, the notation differs a bit from that in |99) . 

4.1. The AKNS Hierarchy. Let L = Jd/dx + Q(x), where 
(4.1.1) J=(* V) ^Q(x)=f ° 



— il \ip(x) 

Note that J 2 = —I, the 2x2 identity matrix, and that JQ + QJ — 0. As mentioned 
in Section 2.3 one may consider the Lax pair (P n+ \, L), where P„+i is a 2x 2-matrix- 
valued differential expression of order n + 1 such that [P rl +i,L] is an operator of 
multiplication. Writing 

n+l 

P n+ i = ^2C n+ i- e (x)L l , 

1=0 

and utilizing that the commutator [Cj,L] may be written as Dj(x)L + Ej(x), the 
condition that [P n+ i,L] is an operator of multiplication yields Dj + \(x) = —Ej(x) 
for j = 0, n, D = 0, and [P n+ i,L] = E n+ i{x). Note that Cj(x) can be expressed 
as 

Cj(x) = kj{x)I + Vj(x)J + Wj(x), j = 0,. .. ,n+ 1, 
where the fcj (x) and (cc) are scalar- valued and the matrices Wj (x) have vanishing 
diagonal elements. Then (' = d/dx) 

D 3 = 2Wj and Ej = ~W Q - QW 3 + v'jl - k'jj + 2v 3 JQ - JW-. 

Since W 3 Q + QW 3 is a multiple of the identity matrix /, D 3+ i + Ej =0 shows that 
k' 3 = for all j = 0, . . . , n + 1, and that the following recursion relation holds, 

(4.1.2) Wo = 0, 

v'jl = W 3 Q + QWj, W j+ x = )-J{W l j-2vjQ), j = 0,...,n+l. 



One concludes 
Next, let 



[P n+1 ,L]=2v n+1 JQ- JW' 



n+l 



K n+ i(z) = ^2 k n +i-jZ j , 



3=0 
n+l 



V n+ i{z,x) = } j v n+ i-j(x)z J 



3=0 
n+l 

y3 



W„+l(z, x) = ^ W n +l-j{x)z j 
3 = 1 

so that the recursion relation ( 4.1.2| ) becomes 

V^ +1 {z,x)I = W n +i(z,x)Q{x) +Q(x)W„ +1 (z,x), 
W n+1 (z,x) = 2V n+1 (z,x)Q(x) -2zJW n+1 (z,x) + J[P n+1 ,L]. 



ELLIPTIC SOLUTIONS OF SOLITON HIERARCHIES 



2!) 



Hence, 

(4.1.3) [P n +x,L] = 2zW n+1 (z,x) + 2V n+1 (z,x)JQ(x) - JW' n+1 {z,x). 

If P n +i is of order n + 1 we define the n th order AKNS equations by 

AKNS„(Q) =Q t - [Pn+uL] =0. 
The first few of these equations are 

Qt = -v Q Q' + 2 Cl JQ, 

Qt = -|j(g" - 2Q 3 ) - cxQ 1 + 2c 2 JQ, 

Qt = - 6Q 2 g') - ciJ(Q" - 2Q 3 ) - c 2 Q' + 2c 3 JQ, 

etc., 

where Vg and Ci,C2,... are arbitrary integration constants. Upon rescaling the t 
variable one may choose Vq = 1. In terms of the AKNS pair (p, g) the homogeneous 
versions (that is, ci = c 2 = • • • = 0) of these equations read 

qt J \qx 

Pt\ _ ^0 f Pxx - 2p 2 q 

q t ) 2 \-q X x + 2pq 2 
Pt\ _ivo_{ Pxxx - Gpqpx 

qt) 4 \-q X xx + 6pqq a 
etc. 

We also mention an interesting scale invariance of the AKNS equations. Suppose 
Q satisfies one of the AKNS equations, that is, AKNS„(Q) = 0. Suppose a^O and 



.4: 



a 

a' 1 



Then Q = AQA also satisfies AKNS„(Q) = 0. We omit the straightforward proof 
which can be found, for instance, in J88| . 

In the particular case of the nonlinear Schrodingcr (NS) hierarchy, where p(x, t) = 
±q(x,t) the matrix A is unimodular, that is, \a\ = 1. 

Note that the KdV hierarchy as well as the modified Korteweg-de Vries (mKdV) 
hierarchy are contained in the AKNS hierarchy. In fact, setting all integration con- 
stants C21+1 equal to zero, the n th KdV equation is obtained from the (2n) th AKNS 
system by the constraint p(x, t) — 1, while the n th mKdV equation is obtained from 
the (2n) th AKNS system by the constraint p(x,t) = ±q(x,t). 

Just as in the KdV case one makes the following observation: Suppose a poly- 
nomial V n+ i whose coefficients are scalar functions and a polynomial W„+i whose 
coefficients arc 2x2 matrix-valued functions with zero diagonal to be given. Fur- 
thermore, assume that V^ +1 (z,x)I — W n +i{z,x)Q(x) + Q(x)W n +i(z, x) and that 
2zW„ + i(z, x) + 2V n+ i(z, x)JQ(x) — JW' n+1 (z, x) is independent of z. Then the co- 
efficients of V n +i (z, x) and W n +i(z, x) define a differential expression P n +i which 



satisfies (|4. 1 .3|) 



Next, suppose that L] — and, without loss of generality, vq ^ 0. Then, 

(4.1.4) (P n+1 -K n+1 {L)f = (JV n+1 (L 1 x) + W n+1 (L,x)) 2 
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= W n+1 (L,x) 2 -V n+1 (L,x) 2 I 

2n+2 

= "ML™, 
m=0 

where 

(x 

)I = ^ (W n+ i-e(x)W n+ i-k(x) — v n+1 -e(x)v n+ i-k{x)I) 

£-\-k—rn 

is a multiple of the identity matrix. Moreover, differentiating (JV„ + i(z, x) + 
W n +i(z, x)) 2 with respect to x yields 

~2V n+1 (z,x)V' l+1 (z,x)I + W„+i(z,x)W' n+1 (z,x) + W' n+1 (z,x)W„+i(z,x) = 0, 

using 

V' l+1 (z,x)I = W„+i(z,x)Q(x) + Q(x)W n+ i(z,x), 
W' n+1 (z,x) = 2V n+ i(z,x)Q - 2zJW n+ i(z,x), and 
= JW n+1 (z,x) + W n+1 {z,x)J. 



Hence the coefficients a m (x) in (4.1.4) may be interpreted as constant scalars. 
Since a 2n +2 = —Vq one infers {P n +i - K n+ i(L)) 2 + R 2n +2(L) = 0, where R 2n +2 is a 
polynomial of degree 2n + 2 with complex coefficients. Hence, if [P n +i , L] — 0, then 
the pair (P n+ i, L) is associated with a hyperelliptic curve of (arithmetic) genus n. 

Next assume that F n : C 2 — > C m is a polynomial of degree n in its first variable 
with scalar meromorphic coefficients. Denote the leading coefficient by —iq(x) and 
let p{x) be another nonzero meromorphic function. Defining 

V n+1 (z,x) = —^—(F'^z.x) + 2izF n (z,x)) 
2q(x) 

and 

q(x)F n (z, x) 



W n+ i(z,x) 



<1( 



(x)\V£ +1 (z,x)+p(x)F n {z,x) 



this implies V^ +1 I = W n +\Q + QW„+i with Q given as in ( 4.1.1 ). Moreover, 
V 2 +1 I — W 2 +1 = i?2n+2(z, x)I, where the scalar i?2n+2(z, x) is given by 

(4.1.5) 

R2n+2(z, x) =_L_ (i^( 2 , xf - 2F n {z, x)F^{z, x) + 4{p{x)q(x) - z 2 )F n {z, x) 2 ) 



+ {2F n (z, x)F' n {z, x) + AizF n {z, x) 2 ) . 



4q(x) 

If this is constant then differentiation with respect to x yields 

W n+1 (z,x)(W^ +1 {z,x) - 2V n+1 (z,x)Q(x)) 

+ {W n+1 {z,x) - 2V n+1 (z,x)Q(x))W n+ x(z,x) = 0. 

Since W' n+1 (z, x) — 2V n +i(z, x)Q(x) has zero diagonal elements, this may be con- 
sidered, for each fixed x, a linear homogeneous equation for the off-diagonal ele- 
ments of W' n+1 {z, x) — 2V n +\(z 1 x)Q(x). This equation has a one-dimensional space 
of solutions, in fact, W! n+1 (z,x) — 2V n +i(z, x)Q(x) — r(x)JW n +i{z,x). Compar- 
ing the leading coefficients yields r(x) = —2z and hence we have shown that 
W' n+1 (z,x) = — 2z JW„ + i(z, x) + 2V n +i(z, x)Q(x). Summarizing, the following the- 
orem holds. 
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Theorem 4.1. If F n : C 2 — > C M is a polynomial of degree n in the first variable 



and meromorphic in the second, and if the expression i?2n+2(z, x) in (4.1.5) is 
independent of X, then there exists a 2 x 2 matrix-valued differential expression 
Pn+i of order n + 1 with leading coefficient J n+2 such that [P n +i, L] = 0. 

We are now ready to define the term algebro-geometric in the context of the 
AKNS hierarchy. 

Definition 4.2. Suppose Q is a 2 x 2 matrix-valued meromorphic function and 
let L be the differential expression L = Jd/dx + Q. Then Q is called an algebro- 
geometric AKNS potential (or simply algebro-geometric) if Q is a solution of some 
equation of the stationary AKNS hierarchy. 

Equivalently, Q is algebro-geometric if any one of the following two conditions 
is satisfied. 

(1) There exists a differential expression P n +i of order n + 1 with leading coeffi- 
cient J n+2 such that [P n+1 ,L] = 0. 

(2) There exists a function F n : C 2 — > C M which is a polynomial in the first 
variable, meromorphic in the second, such that the expression R2 n +2{z, x) in 



(4.1.5) does not depend on x. 



4.2. Periodic AKNS Potentials. Throughout this section we will assume that 
p,q G C 2 (R) are complex- valued periodic functions of period fl > 0. Matrix- valued 
solutions Y of the initial value problem (LY)(x) = zY(x) and Y{x$) = Yq are 
denoted by •, Xq, Yq), which is also the unique solution of the integral equation 

(4.2.1) Y(x) = exp(z(x - x )J)Y + f e zix ~ x ' )J JQ(x')Y(x')dx' . 

Jxo 



In contrast to the Sturm-Liouville case, the Volterra integral equation (4.2.1) 



is not suitable to determine the asymptotic behavior of solutions as z tends to 



infinity. To circumvent this difficulty one can follow Marchenko's approach in [ 153 1 , 
Section 1.4, to obtain the asymptotic expansion 

Hz,x + n,x Q ,i) = ^ e ^)+^- z [ Mxo)MzQ) _ 0e i Z n 

(4.2.2) +0(e |Im(z)|f V 2 ), 
where 

0= / p(t)q(t)dt, 

J Xq 

provided p, q £ C 2 (R). From this result we infer that the entries of <&(•, Xq+Q, Xq, I), 
which are entire functions, have order of growth equal to one whenever q(xo) and 
p(xo) are nonzero. 

As in the scalar case wc introduce 7n(z), the restriction of the translation oper- 
ator Y i— > Y(- + ft) to the two-dimensional space of solutions y{z) of (LY)(x) — 
zY(x). The Floquet multipliers arc then determined as solutions of 

p 2 - ptv(T n (z)) + 1 = 0. 

They are degenerate if and only if p(z) 2 — 1. The points z where this happens 
are the (anti-)periodic eigenvalues which we denote by Ej, j S Z. Their algebraic 
multiplicities are given by ordp. (Tq — 4) and are denoted by p{Ej). The asymptotic 



32 



FRITZ GESZTESY AND RUDI WEIKARD 



behavior of the (anti-)periodic eigenvalues may be determined from (4.2.2). It is 
given by 

E 2j ,E 2j . 1 = ^+0(\j\- 1 ), 

if the eigenvalues are labeled in accordance with their algebraic multiplicity. 
The conditional stability set S(L) of L is again 

S(L) = {zeC\-2< tr(T n («)) < 2}. 

The spectrum of the maximal operator in L 2 (K) 2 associated with L coincides with 
the conditional stability set S(L) of L. One can also prove that the conditional sta- 
bility set S(L) consists of a countable number of regular analytic arcs, the spectral 
bands. At most two spectral bands extend to infinity and at most finitely many 
spectral bands are closed arcs. The point z is a band edge, that is, an endpoint of a 
spectral band, if and only if (tr(7n(z))) 2 — 4 has a zero of odd order. For additional 



results on nonself-adjoint periodic Dirac-type operators, see Tkachenko [214 , [El 
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The boundary value problem defined by the differential expression L and the 
requirement that the first component of a solution vanishes at Xq and xq + £1 will 
(somewhat artificially but in close analogy to the KdV case) be called the Dirichlet 
problem for the interval [xq, xq + £1]. The Dirichlet eigenvalues and their algebraic 
multiplicities are given as the zeros and their multiplicities of the function 

g(z,x ) = (l,0)$(z,x o + fi,j; o ,/)(0,l) t , 

that is, the entry in the upper right corner of <&(z,xq + Cl,Xo,I). The algebraic 
multiplictiy of z as a Dirichlet eigenvalue n{x) is denoted by d(z, x). The quantities 

di(z) = min{c?(z, x) \ x £ R} and d m (z, x) = d{z, x) — di(z) 

will be called the immovable part and the movable part of the algebraic multi- 
plicity d{z, x), respectively. The sum X^zec d m (z, x), which is independent of x, is 
called the number of movable Dirichlet eigenvalues. Asymptotically, the Dirichlet 
eigenvalues are distributed according to 

^(xo) = J -£ + o(\j\- 1 ) 1 jez. 

If q(x) 7^ 0, the function g{-, x) is an entire function with order of growth equal to 
one. Hadamard's factorization theorem then implies g{z,x) = F(z, x)D(z), where 
F comprises the factors depending on x, while D contains the factors independent 
of x. Of course, ord z (_F(-, x)) = d m (z,x) and ord z (-D) = di(z). 

We now turn to the ^-dependence of the function g. The following is the analog 
of equation (3.1.2) and is obtained by a straightforward computation. The function 
F(z, •) satisfies the differential equation 

(4.2.3) q(x){F'(z,x) 2 - 2F{z,x)F" {z,x) + A{p{x)q{x) - z 2 )F{z,x) 2 ) 

+ q'(x)(2F(z, x)F'{z, x) + 4izF{z, xf) = q(x) 3 ({tr(T n (z))f - A)/D{z) 2 . 

Just as in the KdV case we get the following theorem. 

Theorem 4.3. There exists an entire function R such that (tr(l}i(z))) 2 — 4 = R{z) 
D(z) 2 . In particular, p(z) — 2di{z) > for every z£C. 

We define the Floquet deficiency of Q (or L) as def(Q) = deg(i?) e No U {oo}. 
Equation (4.2.3), Theorem 4.1, and the asymptotic distribution of (anti-)periodic 
and Dirichlet eigenvalues allow one to prove the following central theorem. 
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Theorem 4.4. The following statements are equivalent: 

(1) def(Q) = 2n + 2. 

(2) There are n movable Dirichlet eigenvalues. 

(3) There exists a 2 x 2 matrix-valued differential expression P n +i of order n + 1 
and leading coefficient J n+2 which commutes with L. The number n + 1 is the 
smallest integer with that property. The differential expression P„+i satisfies 
P 2 +1 = R2n+2{L) = n ze c{L — z) p ( z ^~ 2di ( z } and hence Q is algebro-geometric. 

We remark that in (3) i?2n+2 is a constant multiple of R. 

4.3. Necessary Conditions for an AKNS Potential to be Algebro-Geome- 
tric. In this section we will prove the following result. 

Theorem 4.5. Suppose Q is an algebro-geometric AKNS potential. Then LY = 
JY' + QY = zy has a meromorphic fundamental system of solutions with respect 
to x for all values of the spectral parameter z G C. 



A proof of this theorem along the lines of the proof of Theorem 3.10 fails since 
one can only show that the exponents of a singularity are half integers rather than 
integers. However, such reasoning can be used to prove the following result. 

Theorem 4.6. Suppose Q is a meromorphic potential coefficient of L. Then the 
equation (LY)(x) = zY(x) has a fundamental system of solutions meromorphic 
with respect to x for all values of the spectral parameter z£C whenever this holds 
for infinitely many values of z. 



The rest of this section is devoted to a proof of Theorem 4.5. Suppose Q is an 
algebro-geometric AKNS potential. According to the characterizations following 
Definition 4.2 there is a function F n : C 2 — > C m which is a polynomial of degree n in 
the first variable, meromorphic in the second, such that the expression i?2n+2(*z, x) 
in (4.1.5) does not depend on x. With the aid of F„ we define (as above) 

V n+1 (z,x) = — n -(F^(z,x) + 2izF n (z,x)), 
2q(x) 

VV n+1 {z,x) - ^ [y; +1 (z,x) +p(x)F n (z,x) 
and 

Pn+l = JV n+1 (L,x)+W n +l(L,x). 

The pair (P n+ i, L) is associated with the hyperelliptic curve 

JC n = {(z, w) :w 2 + R 2n +2{z) = 0} 
of arithmetic genus n, where 

2n+2 

R2n+2(z)I = V n+1 {z,x) 2 I-W n+1 {z,x) 2 = ]J {z - E rn )I. 

m=l 

We emphasize that the points E m are not necessarily distinct. 

If ^ = (ipi^ip^y is a common solution of L^> = z^ and P n +i^> = then, 
considering the first component of P n+ i^, we find 

wipi = (P n +i*)i = [(JV n+1 (z,x) +W n+l (2,i))*] 1 = iV n+1 ipx + W«+i,i,2^2 
= (iV n+ i+W n+hh 20)^i, 
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defining <f> — ^2/^1 ■ We now revert this process and define the meromorphic 
function <fi on JC n x C by 

,,, ^ \ w-iV n+ i(z, x) Wn+1,2,1 
cp{{z,wj,x) = — — - — = — -. 

W n+ i A ,2{z,x) w + iV n+ i{z,x) 

We remark that <f> can be extended to a meromorphic function on the compactifi- 
cation (projective closure) of the affine curve K, n - This compactification is obtained 
by joining two points (the points at infinity) to K, n . 
Next we define 

(4.3.1) V , i(( z ! w )> x i x o) — CX P ( / {—iz + q(x')4>((z,w),x'))dx 

(4.3.2) i/j 2 ((z,w),x,x ) = (j>((z,w),x)ip 1 ((z,w),x,x ) 



where the simple Jordan arc from xq to x in (4.3.1) avoids poles of q and </>. One 
verifies with the help of V 7 [ +1 I = W n+1 Q+QW n+1 and W n+1 = 2V n+ iQ-2zJW n+1 
that 

4>'{{z, w),x) = p(x) — q(x)(j>((z, w),x) 2 + 2iz<fi((z, w),x). 
^From this one obtains that 

•«-).«.*)- ftg:;:!;:;:!; 

is a common solution of L 1 ^ = z^ and P n +\*5! = ui'J'. 

Let m(xo), Hn{%o) denote the zeros of F n (-,xo). One then observes that 
the two branches ^f±(z, ■, xo) — ty((z,±w),-,xo) of the function *$>((z,w),-,xo) 
represent a fundamental system of solutions of Ly = zy for all complex numbers z 
different from E%, E 2n +2, A*i(^o), • ••> Mn(^o), since 

(4.3.3) W(*-(z,x,xo),*+{z,x,x )) - =tj -. r = , v 

VVn+1, 1,2{Z, X ) F n (Z,X ) 

where W(f,g) denotes the determinant of the two columns / and g. 

In the special case where K n is nonsingular, that is, when the points E m are 
pairwise distinct, the explicit representation of ^{{z, w), x, xq) in terms of the Rie- 
mann theta function associated with IC n immediately proves that ty±(z, x, xq) are 
meromorphic with respect to x € C for all z £ C\{E±, E2 n +2, fJ-i{xo), (J-n( x o)}- 
A detailed account of this theta function representation can be found, for instance, 
in Theorem 3.5 of |8S[ ]. In the following we demonstrate how to use gauge trans- 
formations to reduce the case of singular curves JC n to nonsingular ones. 

Let Q be meromorphic on C and introduce A(z,x) = (Q(x) + zI)J which turns 
LY = zY into Y' + AY — 0. Then consider the gauge transformation 

&(z,x) = T(z,x)^(z,x). 

If 

A(z, x) — T(z, x)A(z, x)T(z, x)^ 1 — F'(z, x)T(z, x) -1 , 

then A(z, x) = (Q(x) + zI)J, where Q has zero diagonal elements. Moreover, 
+ A$ = 0, that is, {Jd/dx + Q)# = z^> . Next we make the following ex- 



plicit choice suggested by Konopelchenko [128]. Let z <E C be fixed, ^°'(z,-) — 



(ipi (z, •), ^2 '))* be any solution of Ly = zy, and introduce 
<f>W(z,x)=4°\z,x)/4 0) (z i x). 
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Then one defines 



for z £ C\{z}. 

The upper right entry G\^{z, x, x') of the Green's matrix of L is given by 
n , t\ . ip+,i(z,x,x )ip- A (z,x',x ) , 

G 12 {Z,X,X) = l — — rr, X>X. 

W(V-(Z, ;X0),V+{Z, ;Xq)) 

We want to evaluate Gi t 2(z,x, x') on its diago nal (i .e., where x = x'). Since 
ip+ t i(z, xq, xo)ip- t i (z, Xq, x ) = 1 we obtain from ( 4.3.1 ) and ( 4.3.3 ) 

„ , \ .Wn+i,ia{z,%) -F n (z,x) 

Gi, 2 {z,x,x) = % - = . 

2w 2w 

Next note that 

•), *+(*, •)) = det(r(z, x))W(*-{z, ■),*+(*, ■)) 
and that det(T(z, x)) = —i(z — z)/2. With the help of this fact and some compu- 
tations one finds that the upper right entry G\ 2 {z,x,x) of the Green's matrix of 
Jd/dx + Q is 

V> +) i(z,a;)^L -F n+1 (z,x) 



Gi^(z, x,x) = % 



where F n+ i(-,x) is a polynomial of degree n + 1 with leading coefficient —iq(x). 
Moreover, F n+ i(z,x) satisfies 

q{2F n+1 F'; +l - F l2 +l + 4(z 2 - pq)F 2 +1 ) - q'{2F n+1 F: i+l + 4izF 2 +1 ) 

(4.3.4) =-4g 3 (z-z) 2 i? 2 „ +2 (z). 

Hence Q is an algebro-geometric AKNS potential. 

Now suppose that (z,0) is a singular point of /C„, that is, that z is a zero of 
-R2n+2 of order r > 2. Choose 

^{z,x) = 

V V ' ^1,1,3(2,4) 

Then one may show that F n+ \{-, x) has a zero of order at least 2 at z, that is, 

F n+1 (z,x) = (z - z) s Fn(z,x) 



for some s > 2 and some polynomial jft(-, a;) of degree n = n + 1 — s. From (4.3.4) 
one obtains 

q(2F A F^ ~F l2 + 4(z 2 - pq)F?) - q'{2F h F' fl + AizF 2 ) = -4g 3 i? 2 „_ 2s+4 (z), 
where 

R2n-2s+i{z) =(Z~ Z) 2 ~ 2S R 2n +2{z) 

is a polynomial in z of degree 2n — 2s + 4 e (0, 2n + 2). This proves that Q is 
associated with the curve 

Kk = {(z,w) :w 2 + {z- ~zf- 2s R 2n+ 2{z) = 0}. 

Our choice of 4>^ (z, x) led to a curve /C« which is less singular at (z, 0) than 
K. n , without changing the local structure of the curve elsewhere. By iterating 
this procedure one ends up with a curve which is nonsingular at (z,0). Repeating 
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the procedure for each singular point of IC n then results in a nonsingular curve 
and a corresponding Baker-Akhiezer function ty((z, w),x, Xo) meromorphic with 
respect to x S C using its standard theta function representation (see, e.g., |@). 
Since <^°) — — iV n +\ /Wn+1,1,2 is meromorphic, the gauge transformations and their 
inverses map meromorphic functions to meromorphic functions. Combining these 
results proves Theorem |4.5| . 

A systematic account of the effect of Darboux-type transformations on hyper- 
elliptic curves in connection with the KdV, AKNS, and Toda hierarchies will be 
presented in [j36| . 



4.4. Sufficient Conditions for an AKNS Potential to be Algebro-Geome- 
tric A Characterization of Elliptic Stationary AKNS Solutions. 

Definition 4.7. Let p, q be elliptic functions with a common period lattice. Then 

Q = ( 9 p ^J 9 J is called a Picard-AKNS potential if the equation Jtp' (x)+Q(x)^p(x) — 

ztp(x) has a meromorphic fundamental system of solutions with respect to x for all 
values of the spectral parameter z£C. 



We note that according to Theorem 4.6 it is sufficient to show the existence of a 



meromorphic fundamental system for infinitely many values of z in order to prove 
that Q is a Picard-AKNS potential. 

Just as in the KdV case the following theorem characterizes all elliptic algebro- 
geometric AKNS potentials. 

Theorem 4.8. Let Q — (^i p q 9 ^ withp,q elliptic Junctions with a common period 
lattice. Then Q is an elliptic algebro- geometric AKNS potential if and only if it is 
a Picard-AKNS potential. 



Sketch of proof. The necessity of the criterion is the content of Theorem 4.5. The 



sufficiency follows from Picard's Theorem 2.7 and Theorem 4.4 in the same way as 



in the KdV case. □ 



The transformation 

/ e ax+b Q \ 

allows one to prove the following corollary, which slightly extends the class of 
algebro-geometric AKNS pot entia ls Q considered thus far. Such cases have recently 
been considered by Smirnov pOfJ . 



Corollary 4.9. Suppose 



-~iq(x)e~ 2( - ax + b 1 



(4-4.1) Q(x) = [ w{x y {ax+b) o 

where a, b € C and p, q are elliptic functions with a common period lattice. Then 
Q is an algebro-geometric AKNS potential if J^'(x) + Q(x)^(x) — z^{x) has a 
meromorphic fundamental system of solutions with respect to x for all values of the 
spectral parameter z G C. 
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4.5. Examples. With the exception of the studies by Christiansen, Eilbeck, Enol's- 
kii, and Kostov ||o| and Smirnov [ 200 1 , 1 202 ] , [ 204 ] , not too many examples of elliptic 
solutions (p, q) of the AKNS hierarchy associated with higher (arithmetic) genus 
curves of the type w 2 + R.2n+2{z) = have been worked out in detail. The genus 
n = 1 case has been considered, for example, by Its [121] and Pavlov 177]. More- 
over, examples for low genus n for special cases such as the nonlinear Schrodinger 
and mKdV equation are considered, for instance, in B, p5| , [148 1 , |162|, [156], 
[201 . Examples related to equations of the sine-Gordon- type are discussed in 
[205 , [ 206 1 . The following examples in (4.5.1) - (4.5.4) are algebro-geometric 



172 



197 



AKNS potentials as can be proved using the Frobenius method. For details on this 
procedure see p9~|. 



(4.5.1) 


Q(x) 


(4.5.2) 


Q(x) 


(4.5.3) 


Q(x) 


(4.5.4) 


Q(x) 



= in(((x) 



C(X - W 2 ) - C(^2)) 

(n + l)p(x 











ip'{x 



-3»V(x)/(2ei) 




ip(x 



neN, 

■W2)/(2ei) 


- u> 2 )/ej 




e 2 = 0, 



e 2 = 0. 



Incidentally, if p = 1, as in Example ( [4.5.2] ) then J '5' + Qty = zty is equivalent 
to the scalar equation tp'^ - 



#2 



-z 2 ip2, where ^> = (ipi,^)* and ipi = ip' 2 — izip2 



Therefore, if — q i s an elliptic algebro-geometric potential of the KdV hierarchy then, 
by Theorem 3.12 , ip2 is meromorphic for all values of z. Hence "J is meromorphic for 
all values of z and therefore Q is a Picard-AKNS and hence an algebro-geometric 
AKNS potential. Conversely, if Q is an algebro-geometric AKNS potential with 
p = 1 then — q is an algebro-geometric potential of the KdV hierarchy. In particular, 
q(x) = n(n + l)p(x) is again the class of Lame potentials associated with the KdV 



hierarchy and hence a special case of the material discussed in Section 3.2 
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